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Abstract 



In this paper, we give an elaborate and understandable review of traffic cellular automata 
(TCA) models, which are a class of computationally efficient microscopic traffic flow mod- 
els. TCA models arise from the physics discipline of statistical mechanics, having the goal 
of reproducing the correct macroscopic behaviour based on a minimal description of mi- 
croscopic interactions. After giving an overview of cellular automata (CA) models, their 
background and physical setup, we introduce the mathematical notations, show how to per- 
form measurements on a TCA model's lattice of cells, as well as how to convert these 
quantities into real-world units and vice versa. The majority of this paper then relays an 
extensive account of the behavioural aspects of several TCA models encountered in lit- 
erature. Already, several reviews of TCA models exist, but none of them consider all the 
models exclusively from the behavioural point of view. In this respect, our overview fills 
this void, as it focusses on the behaviour of the TCA models, by means of time-space and 
phase-space diagrams, and histograms showing the distributions of vehicles' speeds, space, 
and time gaps. In the report, we subsequently give a concise overview of TCA models that 
are employed in a multi-lane setting, and some of the TCA models used to describe city 
traffic as a two-dimensional grid of cells, or as a road network with explicitly modelled 
intersections. The final part of the paper illustrates some of the more common analytical 
approximations to single-cell TCA models. 

Key words: cellular automata, traffic flow modelling, tempo-spatial behaviour, 

phase-space diagrams 

PACS: 02.50.-r,45.70.Vn,89.40.-a 



In the field of traffic flow modelling, microscopic traffic simulation has always 
been regarded as a time consuming, complex process involving detailed models 
that describe the behaviour of individual vehicles. Approximately a decade ago, 
however, new microscopic models were being developed, based on the cellular au- 
tomata programming paradigm from statistical physics. The main advantage was 
an efficient and fast performance when used in computer simulations, due to their 
rather low accuracy on a microscopic scale. These so-called traffic cellular au- 
tomata (TCA) are dynamical systems that are discrete in nature, in the sense that 
time advances with discrete steps and space is coarse-grained (e.g., the road is 
discretised into cells of 7.5 metres wide, each cell being empty or containing a 
vehicle). This coarse-graininess is fundamentally different from the usual micro- 
scopic models, which adopt a semi-continuous space, formed by the usage of IEEE 
floating-point numbers [1]. True to the spirit of statistical mechanics, all the TCA 
models discussed in this report do not have a realistic microscopic description of 
traffic flows as their primary intent, but are rather aimed at obtaining a correct 
macroscopic behaviour through their crude microscopic description. Such an ap- 
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proach would involve more human- oriented aspects such as those found in socio- 
economic, behavioural, and psychological sciences. Due to large lack of knowledge 
about the manner in which human beings operate in a traffic system, traffic engi- 
neers currently stick with this higher-level scientific approach. As such, they are 
able to positively capture the first- and second-order macroscopic effects of traf- 
fic streams. TCA models are very flexible and powerful, in that they are also able 
to capture all previously mentioned basic phenomena that occur in traffic flows 
[2,3]. In a larger setting, these models describe self-driven, many-particle systems, 
operating far from equilibrium. And in contrast to strictly gaseous analogies, the 
particles in these systems are intelligent and able to learn from past experience, 
thereby opening the door to the incorporation of behavioural and psychological 
aspects [4,5,6]. 

The cellular automata approach proved to be quite useful, not only in the field of ve- 
hicular traffic flow modelling, but also in other fields such as pedestrian behaviour, 
escape and panic dynamics, the spreading of forest fires, population growth and 
migrations, cloud formation, material properties (corrosion, cracks, creases, peel- 
ing et cetera), ant colonies and pheromone trails, . . . [7,8,9,10,11]. It is now fea- 
sible to simulate large systems containing many 'intelligent particles', such that 
is it possible to observe their interactions, collective behaviour, self-organisation, 
...[12,13,7,6,14,15,16,17] 

In this report, we provide a detailed description of the methodology of cellular 
automata applied to traffic flows. We first discuss their background and physical 
setup, followed by an account of the mathematical notations we adopt. The remain- 
ing majority of this report extensively discusses the behavioural aspects of several 
state-of-the-art TCA models encountered in literature (our overview distinguishes 
between single-cell and multi-cell models). The report concludes with a concise 
overview of TCA models in a multi-lane setting, and TCA models used to describe 
two-dimensional traffic (e.g., a grid for city traffic). We end with a description of 
several common analytical approximations to single-cell TCA models. 

I Note that aside from our phenomenological discussion of different TCA 
models, we refer the reader to the work of Chowdhury et al. [3], Santen [18], 
and Knospe et al. [19] for more theoretically- and quantitatively-oriented 
overviews. 



1 Background and physical setup for road traffic 

In this section, we give a brief overview of the historic origins of cellular automata, 
as they were conceived around 1950. We subsequently describe which main ingre- 
dients constitute a cellular automaton: the physical environment, the cells' states, 
their neighbourhoods, and finally a local transition rule. We then move on to a 
general description on how cellular automata are applied to vehicular road traffic, 
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discussing their physical environment and the accompanying rule set that describes 
the vehicles' physical propagation. 

1 . 1 Historic origins of cellular automata 

The mathematical concepts of cellular automata (CA) models can be traced back 
as far as 1948, when Johann Louis von Neumann introduced them to study (living) 
biological systems [20]. Central to von Neumann's work, was the notion of self- 
reproduction and theoretical machines (called kinematons) that could accomplish 
this. As his work progressed, von Neumann started to cooperate with Stanislaw 
Martin Ulam, who introduced him to the concept of cellular spaces. These de- 
scribed the physical structure of a cellular automaton, i.e., a grid of cells which 
can be either 'on' or 'off [21,22]. Interestingly, Alan Mathison Turing proposed 
in 1952 a model that illustrated reaction-diffusion in the context of morphogene- 
sis (e.g., to explain the patterns of spots on giraffes, of stripes on zebras, . . . ). His 
model can be seen as a type of continuous CA, in which the cells have a direct 
analogy with a simplified biological organism [23]. 

In the seventies, CA models found their way to one of the most popular applications 
called 'simulation games', of which John Horton Conway's "Game of Life" [24] 
is probably the most famous. The game found its widespread fame due to Martin 
Gardner who, at that time, devoted a Scientific American column, called "Math- 
ematical Games", to it. Life, as it is called for short, is traditionally 'played' on 
an infinitely large grid of cells. Each cell can either be 'alive' or 'dead'. The game 
evolves by considering a cell's all surrounding neighbours, deciding whether or not 
the cell should live or die, leading to phenomenon called 'birth', 'survival', and 
'overcrowding' (or 'loneliness'). An example of a Life game board can be seen in 
Fig. 1. Typical of Life, is the spawning of a whole plethora of patterns or shapes, 
having illustrious names such as gliders, guns, space ships, puffers, beehives, oscil- 
lators, . . . The Game of Life is now all about how these shapes evolve, and whether 
or not they die out or live indefinitely (either by remaining stationary or moving 
around). 

The widespread popularisation of CA models was achieved in the eighties through 
the work of Stephen Wolfram. Based on empirical experiments using computers, 
he gave an extensive classification of CA models as mathematical models for self- 
organising statistical systems [21,26]. Wolfram's work culminated in his mammoth 
monograph, called A New Kind of Science [26]. In this book, Wolfram related cel- 
lular automata to all disciplines of science (e.g., sociology, biology, physics, math- 
ematics, . . . ). Despite the broad range of science areas touched upon, Wolfram's 
book has received its share of criticism. As an example of this, we mention the 
comments of Gray, who points out that Wolfram's results suffer from a rigourous 
mathematical test. As a consequence, the physical examples in his book are deemed 
either uncheckable or unconvincing. Gray's final critique is that ". . . he [Wolfram] 
has helped to popularise a relatively little-known mathematical area (CA theory), 
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Fig. 1. An example of the Game of Life, with a rectangular grid of cells. Live cells are 
coloured black, whereas dead cells remain white. The image shows a snapshot during the 
game's course, illustrating many different shapes to either die out, or live indefinitely by 
remaining stationary or moving around (image adapted from [25]). 

and he has unwittingly provided several highly instructive examples of the pitfalls 
of trying to dispense with mathematical rigour" [27]. However, with respect to their 
computational power, CA models can emulate universal Turing machines within the 
theories of computation and complexity. Recently, Chua took Wolfram's empirical 
observations one step further, proving that some of the CA models are capable of 
Turing universal computations. He furthermore introduced the paradigm of cellu- 
lar neural networks (CNN), which provide a very efficient method for performing 
massive parallel computations, and are a generalisation of cellular automata [28]. 

Finally, an important step in this direction, is Bill Gosper's proof that the Game 
of Life is computationally universal, i.e., it can mimic arbitrary algorithms [29]. 
Notably, one of the most profound testimonies related to this concept, is the work 
of Konrad Zuse and Edward Fredkin at the end of the sixties. Their Zuse-Fredkin 
thesis states that "The Universe is a cellular automaton" , and is based on the as- 
sumption that the Universe's physical laws are discrete in nature [30,31,32]. This 
latter statement was also conveyed by Wolfram in his famous CA compendium 
[26]. 

1.2 Ingredients of a cellular automaton 

From a theoretical point of view, four main ingredients play an important role in 
cellular automata models [33,22,34]: 

(1) The physical environment 

This defines the universe on which the CA is computed. This underlying structure 
consists of a discrete lattice of cells with a rectangular, hexagonal, or other topol- 
ogy (see Fig. 2 for some examples). Typically, these cells are all equal in size; the 
lattice itself can be finite or infinite in size, and its dimensionality can be 1 (a linear 
string of cells called an elementary cellular automaton or EC A), 2 (a grid), or even 
higher dimensional. In most cases, a common — but often neglected — assump- 
tion, is that the CA's lattice is embedded in a Euclidean space. 
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Fig. 2. Some examples of different Euclidean lattice topologies for a cellular automaton in 
two dimensions. Left: rectangular. Middle: triangular/isometric. Right: hexagonal. 

(2) The cells' states 

Each cell can be in a certain state, where typically an integer represents the number 
of distinct states a cell can be in, e.g., a binary state. Note that a cell's state is not re- 
stricted to such an integer domain (e.g., Z2), as a continuous range of values is also 
possible (e.g., R + ), in which case we are dealing with coupled map lattices (CML) 
[35,36]. We call the states of all cells collectively a CA's global configuration. This 
convention asserts that states are local and refer to cells, while a configuration is 
global and refers to the whole lattice. 

(3) The cells' neighbourhoods 

For each cell, we define a neighbourhood that locally determines the evolution of 
the cell. The size of neighbourhood is the same for each cell in the lattice. In the 
simplest case, i.e., a ID lattice, the neighbourhood consists of the cell itself plus its 
adjacent cells. In a 2D rectangular lattice, there are several possibilities, e.g., with 
a radius of 1 there are, besides the cell itself, the four north, east, south, and west 
adjacent cells {von Neumann neighbourhood), or the previous five cells as well as 
the four north-east, south-east, south-west, and north-west diagonal cells {Moore 
neighbourhood); see Fig. 3 for an example of both types of neighbourhoods. Note 
that as the dimensionality of the lattice increases, the number of direct neighbours 
of a cell increases exponentially. 

(4) A local transition rule 

This rule (also called function) acts upon a cell and its direct neighbourhood, such 
that the cell's state changes from one discrete time step to another (i.e., the system's 
iterations). The CA evolves in time and space as the rule is subsequently applied 
to all the cells in parallel. Typically, the same rule is used for all the cells (if the 
converse is true, then the term hybrid CA is used). When there are no stochastic 
components present in this rule, we call the model a deterministic CA, as opposed 
to a stochastic (also called probabilistic) CA. 

As the local transition rule is applied to all the cells in the CA's lattice, the global 
configuration of the CA changes. This is also called the CA's global map, which 
transforms one global configuration into another. This corresponds to the notion of 
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Fig. 3. Two commonly used two dimensional CA neighbourhoods with a radius of 1: the 
von Neumann neighbourhood (left) consisting of the central cell itself plus 4 adjacent cells, 
and the Moore neighbourhood (right) where there are 8 adjacent cells. Note that for one 
dimensional CAs, both types of neighbourhoods are the same. 

computing a function in automata theory, see also section 2.1. Sometimes, the CA's 
evolution can be reversed by computing past states out of future states. By evolving 
the CA backwards in time in this manner, the CA's inverse global map is computed. 
If this is possible, the CA is called reversible, but if there are states for which no 
precursive state exists, these states are called Garden of Eden (GoE) states and the 
CA is said to be irreversible. 

Finally, when the local transition rule is applied to all cells, its global map is com- 
puted. In the context of the theory of dynamical systems, this phenomenon of lo- 
cal simple interactions that lead to a global complex behaviour (i.e., the sponta- 
neous development of order in a system due to internal interactions), is termed 
self-organisation or emergence. 

Whereas the previous paragraphs discussed the classic approach to CA models, the 
following sections will exclusively focus on vehicular traffic flows, leading to traffic 
cellular automata (TCA) models: section 1.3 discusses the physical environment on 
which these TCA models are based, and section 1 .4 deals with their accompanying 
rule set that determines the vehicular motion. 



1.3 Road layout and the physical environment 



When applying the cellular automaton analogy to vehicular road traffic flows, the 
physical environment of the system represents the road on which the vehicles are 
driving. In a classic single-lane setup for traffic cellular automata, this layout con- 
sists of a one-dimensional lattice that is composed of individual cells (our descrip- 
tion here thus focuses on unidirectional, single-lane traffic). Each cell can either be 
empty, or is occupied by exactly one vehicle; we use the term single-cell models to 
describe these systems. Another possibility is to allow a vehicle to span several con- 
secutive cells, resulting in what we call multi-cell models. Because vehicles move 
from one cell to another, TCA models are also called particle-hopping models [37]. 



7 



An example of the tempo- spatial dynamics of such a system is depicted in Fig. 4, 
where two consecutive vehicles i and j are driving on a ID lattice. A typical dis- 
cretisation scheme assumes AT = 1 s and AX = 7.5 m, corresponding to speed 
increments of AV = AX/ AT = 27 km/h. The spatial discretisation corresponds 
to the average length a conventional vehicle occupies in a closely packed jam (and 
as such, its width is neglected), whereas the temporal discretisation is based on a 
typical driver's reaction time and we implicitly assume that a driver does not react 
to events between two consecutive time steps [38]. 
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Fig. 4. Schematic diagram of the operation of a single-lane traffic cellular automaton 
(TCA); here, the time axis is oriented downwards, the space axis extends to the right. The 
TCA's configuration is shown for two consecutive time steps t and t+l, during which two 
vehicles i and j propagate through the lattice. 

With respect to the layout of the system, we can distinguish two main cases: closed 
versus open systems. They correspond to periodic (or cyclic) versus open boundary 
conditions. The former is usually implemented as a closed ring of cells, sometimes 
called the Indianapolis scenario, while the latter considers an open road. This last 
type of system, is also called the bottleneck scenario. The name is derived from 
the fact that this situation can be seen as the outflow from a jam, where vehicles 
are placed at the left boundary whenever there is a vacant spot. Note that, in closed 
systems, the number of vehicles is always conserved, leading to the description of 
number conserving cellular automata (NCCA) [39]. 

1.4 Vehicle movements and the rule set 

The propagation of the individual vehicles in a traffic stream, is described by means 
of a rule set that reflects the car-following and lane-changing behaviour of a traf- 
fic cellular automaton evolving in time and space. The TCA's local transition rule 
actually comprises this set of rules. They are consecutively applied to all vehicles 
in parallel (called a parallel update). So in a classic setup, the system's state is 
changed through synchronous position updates of all the vehicles: for each vehi- 
cle, the new speed is computed, after which its position is updated according to 
this speed and a possible lane-change manoeuvre. Note that there are other ways to 
perform this update procedure, e.g., a random sequential update (see section 3.2.4). 
Because time is discretised in units of AT seconds, an implicit reaction time is 
assumed in TCA models. It is furthermore assumed that a driver does not react to 
events between consecutive time steps. 
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For single-lane traffic, we assume that vehicles act as anisotropic particles, i.e., 
they only respond to frontal stimuli. So typically, the car-following part of a rule 
set only considers the direct frontal neighbourhood of the vehicle to which the rules 
are applied. The radius of this neighbourhood should be taken large enough such 
that vehicles are able to drive collision-free. Typically, this radius is equal to the 
maximum speed a vehicle can achieve, expressed in cells per time step. 

From a microscopic point of view, the process of a vehicle following its prede- 
cessor is typically expressed using a stimulus-response relation [1]. Typically, this 
response is the speed or the acceleration of a vehicle; in TCA models, a vehicle's 
stimulus is mainly composed of its speed and the distance to its leader, with the 
response directly being a new (adjusted) speed of the vehicle. In a strict sense, this 
only leads to the avoidance of accidents. Some models however, incorporate more 
detailed stimuli, such as anticipation terms. These forms of 'anticipation' only take 
leaders' reactions into account, without predicting them. When these effects are 
taken into account together with a safety distance, strong accelerations and abrupt 
braking can be avoided. Hence, as the speed variance is decreased, this results in a 
more stable traffic stream [40,41,42]. 

To conclude this section, we note that a TCA model can also be derived from a 
so-called Gipps car-following model. All speeds in this Gipps model are directly 
computed from one discrete time step to another [1]. If now the spatial dimension 
is also discretised (a procedure called coarse graining), then this will result in a 
TCA model. 



2 Mathematical notation 

In this section, we give an overview of the mathematical notation adopted through- 
out this report. The focus will be on the variables in TCA models, the measurements 
that can be done on a TCA model's lattice, and their conversion to real-world units. 
We first take a look at the notation that is commonly used in automata theory, from 
which cellular automata sprung. 

2.1 Classic notation based on automata theory 

Let us first briefly present the notation for cellular automata models, adopted in 
spirit of automata theory. As mentioned in section 1, a CA model represents a 
discrete dynamic system, consisting of four ingredients: 

CA=(C,E,tf,6), CD 

where the physical environment is represented by the discrete lattice C and the set 
of possible states denoted by E. Each i th cell of the lattice, has at time step t a state 
<7j(i) G E. Furthermore, the associated neighbourhood with this cell is represented 
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by Afi(t), i.e., a (partially) ordered set of cells. Finally, the local transition rule is 
represented as: 

5 : E |AA| — > E : (J (7,(0 i— > ^(t + 1). (2) 

Equation (2) shows that the state of the ith cell at the next time step t + 1 is com- 
puted by 5 based on the states of all the cells in its neighbourhood at the current 
time step t. In the previous equation, \M\ represents the number of cells in this 
neighbourhood, which is taken to be invariant with respect to time and space. Note 
that the local transition rule is commonly given by a rule table, where the output 
state is listed for each possible input configuration of the neighbourhood. Given the 
sizes of E and Af, the total number of possible rules equals: 

where each of the possible configurations of a cell's neighbourhood is mapped 
to the number of possible states a cell can be in. 

Considering the ordered set of all the states of all cells collectively at time step t, a 
CA's global configuration is obtained as: 

C(t) = (J ( 4 ) 
jec 

with C(t) E E £ where the latter refers to the set of all possible global configurations 
a CA can be in (also called its phase space). Sometimes, such a global configuration 
C(t) is also represented by its characteristic polynomial (i.e., generating function) 
[43]: 

3=0 

If we now apply the local transition rule to all the cells in the CA's lattice, the next 
configuration of the CA can be computed by its induced global map: 

G:E £ — >E £ :C(t)i — >C(t + l). (6) 

Note that if the CA is reversible, the inverse global map G' 1 can be computed. As 
the CA evolves in time and space, the global map is iterated from a certain initial 
configuration C(0) at t = 0, leading to the following sequence of configurations: 

C(0) - G(C(0)) - G 2 (C(0)) - G 3 (C(0)) - • • • (7) 
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The above sequence is called the trajectory of the initial configuration C(0) under 
the global map G, and we denote it by: 

T m \ G = {G n (C(0)) | n e N}. (8) 

When this trajectory is periodic or chaotic, we use the terminology forward orbit 
and denote it by Cc(o)|g- Similarly, the backward orbit (i.e., the reverse trajectory) 
is denoted by O c ^ a , , where we specify a certain global configuration C(t) at time 
step t under the inverse global map G~ l . 

2.1.1 Classification of CA rules 

Computing the global map G is rather difficult, as it may require many or even 
an infinite amount of iterations in order to obtain the trajectories. In practice, the 
system's lattice size should be taken infinitely large, but even only considering 1000 
cells of a binary elementary cellular automaton (ECA) would increase the size of 
the search space of global configurations to 2 10()0 10 300 . 

A more intuitive methodology, is to observe a CAs tempo-spatial behaviour, i.e., its 
evolution on the lattice in the course of time. To this end, Stephen Wolfram empir- 
ically studied many configurations of binary ECA rules, with a neighbourhood of 
three cells. According to equation (3), this amounts to 2 = 256 different rules. In 
1984, based on this research, Wolfram conjectured four distinct universality classes 
[44]: 

Class I 

These CA evolve after a finite number of iterations to a unique homogeneous 
state, i.e., a limit point. 
Class II 

These CA generate regular, periodic patterns, i.e., entering a limit cycle. 
Class III 

CAs in this class evolve to aperiodic patterns, independent of the initial config- 
uration; their trajectories in the configuration space lie on a chaotic attractor. 
Class IV 

This class encompasses all the CAs that seem to behave in a complex way, 
with features such as propagating structures, long transients; they are thought 
to have the capability of universal computation. 

Although Wolfram's classification scheme is widely adopted, it still remains a ten- 
tative result as he himself states [26]. Note that the type of classification he provides 
is phenotypic, in the sense that it is based on observed behaviour, whereas a geno- 
typic classification would be based on the intrinsic structure of the rules in each 
class. 

Despite these observations, classification still remains a difficult task as is evi- 
denced by the ongoing research in dynamical systems. Other attempts at classi- 
fication of ECA rules include the following. Firstly, Culik and Yu gave a formali- 
sation of Wolfram's classes [45]. Secondly, Li and Packard studied the structure of 
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the ECA rule space according to a certain distance metric, resulting in five classes 
[46]. Then, Braga et al. identified three classes based on the growth of patterns 
observed in CA models [47]. Next, Wuensche used a whole arsenal of local mea- 
sures to automatically create complex rules, thereby classifying the rule space for 
the CAs' dynamics [48]. Furthermore, Dubacq et al. classified CA models based 
on their algorithmic complexity by measuring the information content of the local 
transition rule [49]. And finally, Fates who used a macroscopic parameter, i.e., the 
density of l's, to separate chaotic ECA rules from non-chaotic ones [50]. 

2.1.2 An example of a CA 

To end this section, let us give some definitions of a one dimensional, infinitely 
large, binary state CA with a neighbourhood of radius 1 : 



C = Z d (withd=l), (9) 

£ = Z 2 = {0,1}, (10) 

M = {*-i,M + i}, (ID 

8(i,t) : Z 3 2 — > Z 2 

: {a i _ 1 (t),<r i (*),(T i+1 (0}'— + (12) 
G(C{t)) : Zf — ► Zf 

: C(t) . — » 1). (13) 



Note that in equation (11), we assume that the i th cell's neighbourhood is repre- 
sented by integer indices (i.e., the cells form a totally ordered set). This alleviates 
the need for an explicit representation of the cells themselves, as it is now sufficient 
to work with the cells' indices and states. The transition rule 5 in equation (12) 
takes as its arguments a cell's index i and current time step t, but operates on the 
states of this cell's neighbourhood. The global map in equation (13) operates on the 
global configuration of the CA at time step t. 

2.2 Basic variables and conventions 

Conforming to the setup and notation discussed in the previous sections, we denote 
a TCA's discrete lattice by C (for the remainder of this section, we assume a rect- 
angular lattice). This lattice physically represents the road on which vehicles will 
drive in a TCA model. It consists of L c lanes, each of which has K c cells, so in 
total there are L c x K c cells in the lattice (L c , K c E N ). Each cell can either be 
empty, or occupied with a single vehicle that spans one or more consecutive cells. 
An example of a lattice containing several vehicles, can be seen in Fig. 5. 

Based on the microscopic vehicle characteristics of a vehicle's space headway, 
space gap, length, time headway, time gap, and occupancy time, we propose to 
use the following set of definitions for multi-lane vehicular road traffic flows that 
are heterogeneous (in the sense of having different vehicle lengths) [51]: 
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Fig. 5. A portion of the lattice C at a certain time step; it has Lc = 3 lanes, containing six 
vehicles. The central vehicle i has a space headway h Si = 6 cells, consisting of a space gap 
g Si = 4 cells and its length k = 2 cells. There are four other space gaps to be considered 
when the neighbouring lanes are taken into account: g l s 'f (left-front), g l s f (left-back), gl'f 
(right-front), and glf (right-back), equalling 6, 4, 2 and 2 cells, respectively. 

g l a >f = x l ; f -Xi-li, (14) 

9 s/ x i x i ^ii (15) 
Lb Lb il,b /i £\ 

9 si = x i~ x i ~h , (16) 

r.b r,b rr,b /m\ 

9s\ =Xi-Xi' -k' , (17) 

for which we assume that a vehicle's position is denoted by the cell that contains 
its rear bumper. For the example in Fig. 5, the left and right frontal and backward 
space gaps of the central vehicle i are 6, 4, 2 and 2 cells, respectively (all these space 
gaps thus represent effective distances, corresponding to the number of empty cells 
between vehicles). Similar definitions hold for the space headways h^f, h r g 'f, h l £, 
and h r s f, i.e., the vehicle lengths in the right-hand sides of equations (14) - (17) are 
dropped. Derivations for the time gaps g\'f , g r t f , g l t f, and g\f , and time headways 
h l t 'f, h r t f, h l t f, and h r t f are analogous. 

Discriminating between frontal an backward neighbours in the adjacent lanes to the 
i th vehicle, is done based on their positions, i.e.: 

{x^\xl' b }< Xl <{xY,xl J }. (18) 

According to equation (18), a vehicle that is driving alongside in an adjacent lane 
to the i th vehicle, will be considered as a backward neighbour as long as its rear 
bumper is located strictly behind the rear bumper of the 2 th vehicle (even if this 
neighbour has a large length that 'sticks out' in front of the i th vehicle). 

Under the above set of assumptions, we can now write the conditions for a success- 
ful lane change (i.e., a possible gap acceptance) as the following constraints: 

9 l s/>0 A ^f>0(left lane change), (19) 
9 r s J > A g r s ' b > (right lane change). (20) 
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With respect to the domains of all variables, we note that all vehicle lengths, space 
gaps, and space headways are expressed as integers, or more specifically: 

k,h Si ,h l *,h r s feN , 
g Si ,h^,h r ;feN, 

In contrast to this, the occupancy times, time headways, and time gaps are not 
restricted to the domain of integers, i.e.: 

Pi ,h ti ,h l t f,h r t ?eR+, 
gti ,h l t 'f,h r t feR + , 

Lb r,b Lf r, f _ jrr, 

g ti , 9t\ , 9t; , 9u e K. 

For example, the occupancy time p { as defined by p { = lijvi [51], corresponds to 
the time a vehicle 'spends' in its own cells. 

To conclude, each vehicle % in the lattice has an associated speed v i E N (expressed 
in cells per time step AT), which is bounded by a maximum speed t> max E N . 
For example, if we set AT = 1.2 s, AX = 7.5 m, and t> max = 5 cells/time step, 
then Vi E {0, . . . , 5} which corresponds to a maximum of 5 x AX/AT = 5 x 
7.5 m/s 4- 1.2 s = 31.25 m/s = 112.5 km/h. As can be seen in this derivation, we 
only consider positive speeds in our models, i.e., vehicles always move forward. 

2.3 Performing macroscopic measurements 

The previously discussed quantities are all microscopic traffic stream characteris- 
tics. In this section, we reconsider the macroscopic quantities densities, flows, and 
mean speeds [51]. As we now have to measure these quantities on a TCAs lattice 
C, we present three possibilities for obtaining the data points: 

• by performing local measurements with an artificial loop detector of finite length 
(open and closed systems), 

• by performing global measurements on the entire lattice (closed system), 

• and by performing local measurements with an artificial loop detector of unit 
length (open and closed systems). 

In the following three sections, we give detailed derivations of each of these mea- 
surement techniques. Locally measured quantities are indicated by a T subscript, 
whereas globally measured ones are indicated by an l g' subscript. A temporal and 
spatial discretisation of respectively AT (in seconds) and AX (in metres) is im- 
plicitly assumed. 
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For all following techniques, we assume an integer measurement period of 
T mp time steps. Thus, aggregating data into intervals of 60 seconds with 
AT = 1.2 s, requires a measurement period of: 



mp 



60 



1.2 



50 time steps. 



(21) 



Furthermore, densities are expressed in vehicles per cell, flows in vehicles 
per time step, and space-mean speeds in cells per time step. 



2.3.1 Local measurements with a detector of finite length 



In this section, we deal with an artificial loop detector of finite length K\ d e No, 
located in a single lane. Note that typically, K\ A > f max , so as to ensure that no 
vehicles can 'skip' the detector between consecutive time steps. The first step in 
our approach for performing these measurements, is based on obtaining local mea- 
surements of the density and flow for such a spatial measurement region at a certain 
time step t [51]. Once these are known, the space-mean speed can be derived using 
the fundamental relation of traffic flow theory q = kv s [51]: 



h(t) 



N(t) 
Km ' 

1 N(t) 



9i(*) = ir 



Id i=1 



4 



Qi(t) 



N(t) 



h(t) N(t) ^ 



(22) 
(23) 

(24) 



where we assumed N(t) vehicles are present at time t in the loop detector's seg- 
ment. The density and flow measurements of consecutive time steps are now tem- 
porally averaged over subsequent spatial measurement regions. In similar fashion 
as before, the space-mean speed is derived using the previously mentioned funda- 
mental relation: 
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i Trap i Trap 

k = ■=- E fc < (*) = Tr—fT- E^(*)> (25) 
ft = TjT- E «(*) = ^-sr- E E (26) 

_ T mp ^(*) / T mp 

«.. = f = EE«i(*) / E*(*), (2 7 ) 

^ t=l i=l / t=l 

T mp , AT(t) / T mp 

Jmp / T m p 

= E^(*) / E^to- (28) 

t=i / t=i 

Our derivations for fcj and as outlined above, also correspond to the generalised 
definitions of density and flow, defined as the total time spent, respectively the total 
distance travelled, divided by the area of the measurement region (which corre- 
sponds to T mp x K\ d ). Furthermore, note that the last equation (28) essentially is a 
weighted mean of the local space-mean speeds v Sl (t) at each time step t, with the 
number of vehicles N(t) as weights. 



2.3.2 Global measurements on the entire lattice 

For the global measurements, we consider N vehicles that are driving in a closed 
single-lane system, i.e., with a length of K c cells (the extension to multi-lane traf- 
fic is straightforward). As a consequence, the global density k g remains constant 
during the entire measurement period. The derivations of the equations for k g , q g , 
and v Sg , are completely equivalent to those of the previous section 2.3.1, but now 
with K\ A = K c : 



N 



Tmp N 



T mp Kc t=l i=1 



EE^)> 



4 



Tmp N 



K 9 J mp ^ t=l i= i 

1 Tmp 1 N 

E^E ««(*), 



(29) 
(30) 

(3D 



T N ^ N 

-i mp JV t=1 iv i=l 

^ Tmp 

^- E ■ 

*=i 



(32) 
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Note that, for single-cell TCA models, the global density computed with equa- 
tion (29) actually corresponds to the macroscopic characteristic called occupancy 
p [51]. For multi-cell models, the number of vehicles is in general less than the 
number of occupied cells. 



2.3.3 Local measurements with a detector of unit length 

The third technique for measuring macroscopic traffic flow characteristics on a 
TCA's lattice, bears perhaps the closest resemblance to reality: it is based on an 
artificial loop detector with unit length, i.e., K ld = 1 cell. The loop detector now 
explicitly counts all the vehicles that pass it at each time step AT during the mea- 
surement period T mp . 

This type of measurement corresponds to a point measurement in a temporal mea- 
surement region. Because of this, the appropriate method for computation is dif- 
ferent from the one used in the previous two sections: we now first compute the 
local flow, and the local space-mean speed, both for single-lane traffic. The local 
density is then derived according to the previously mentioned fundamental relation, 
resulting in the following set of equations: 



(33) 
(34) 




(35) 



in which N now denotes the number of vehicles that have passed the detector during 
the measurement period T mp . Because the detector physically occupies one cell and 
because a vehicle has to 'drive by' in order to get counted, this means that stopped 
vehicles are ignored: only moving vehicles are counted. Note that, as opposed to 
the previous two techniques, the above measurements no longer denote temporal 
averages. And because we are working with a temporal measurement region, we 
have to take the harmonic average of the vehicles' speeds Vi in order to obtain the 
local space-mean speed v Sl [51]. 



2.4 Conversion to real-world units 



Converting between TCA and real-world units seems straightforward, as we only 
need to suitably multiply with or divide by the temporal and spatial discretisations 
AT and AX, respectively. However, problems arise due to the discrete nature of a 
TCA model, involving some intricacies with respect to coordinate systems and their 
associated units. For example, as defined in section 2.2, a vehicle i's space headway 
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h Si is always an integer, expressing the number of cells. The same holds true for 
its space gap g Si and length Zj. The difficulty now lies in the fact that fractions of 
cells are not representable in our definition of a TCA model. Keeping in mind that 
h Si — g Si + k [51], and noting that h s . e N , it follows that g s . + Zj > 0, which 
means that either g Si ^ and/or Zj ^ 0. 

As a solution, we therefore adopt throughout this report the convention that, without 
loss of generality, a vehicle's length Zj > 1 cell (which agrees perfectly with our 
earlier definitions in section 2.2). Consequently, when a vehicle % is residing in a 
compact jam (i.e., 'bumper-to-bumper' traffic), its space headway h Si = I cells 
and its space gap g Si = cells. Our convention thus gives a rigourous justification 
to formulate the TCA's update rules more intuitively using space gaps, because as 
already stated in section 1 .4, the rules in a TCA rule set are typically not expressed 
in terms of space headways, but rather in terms of speeds and space gaps (i.e., the 
distance to the leading vehicle). 

In a similar fashion, time headways, time gaps, and occupancy times represent mul- 
tiples of the temporal discretisation AT. But note that, as explained before in sec- 
tion 2.2, these are however no longer constrained to integer values. 

In the following two sections, we explain how to convert between coordinate sys- 
tems of TCA models and the real world. All common variables (e.g., h Si ) are ex- 
pressed in TCA units, except for their 'primed' counterparts (e.g., h' s .), which are 
expressed in real-world units. The conversions will be done with respect to the 
following conventions: 

• TCA model 

• h Si , g Si , and k are dimensionless integers, denoting a number of cells, 

■ h ti , g ti , and pi are dimensionless real numbers, denoting a fractional multiple 
of a time step, 

• ki and kg are real numbers, expressed in vehicles/cell, 

• qi and q g are real numbers, expressed in vehicles/time step, 

■ and Vi, v sr and v Sg are real numbers, expressed in cells/time step. 

• Real world 

■ AX, h' s ., g' s ., and l[ are real numbers, expressed in metres, 

• AT, h' t ., g' t ., and p\ are real numbers, expressed in seconds, 

• k[ and k' g are real numbers, expressed in vehicles/kilometre, 

• q[ and q' g are real numbers, expressed in vehicles/hour, 

• and v[, v 1 , and v' s are real numbers, expressed in kilometres/hour. 



2.4.1 From a TCA model to the real world 

Under the previously mentioned convention that Zj e N , we can write the conver- 
sions of the microscopic characteristics related to the space and time headways and 
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gaps, and the vehicle lengths and occupancy times, in a straightforward manner: 



h'=h Si -AX, g'=g Si -AX, l' t = h-AX, 



h' ti = h ti ■ AT, g' ti = g ti ■ AT, p\ = Pi ■ AT. 



(36) 



Related to equations (36), there is a small but important detail that is easily over- 
looked: we can not just convert between g Si , g' s ., and /■ without making some 
assumptions. Because we adopted the convention that k > I cell, it follows that 
/• > AX. So it is not possible to take the real length of a vehicle smaller than the 
spatial discretisation, because we assumed that the spatial units of a TCA model 
are all integer values. 



The conversions for the macroscopic traffic stream characteristics densities, flows, 
and space-mean speeds, as well as the microscopic vehicle speed, are as follows: 



k' = k- 



q = q- 



1000 

AX ' 

3600 



v* ■ 3.6 ■ 



AX 
AT 



(37) 



To keep the previous equations clear, we have dropped the subscripts denoting 
global and local measurements. 
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It is interesting to see what happens at the jam density, i.e., the maximum 
density when all cells in the lattice are occupied. As all vehicles are standing 
still bumper-to-bumper, the associated space gap at this density, equals zero. 
Computing the space headway, results in h Si — + U. By virtue of the fact 
that density is inversely proportional to the average space headway [51], 
we can cast this space headway into a density, e.g., for a single-cell TCA 
model: kj — h Sj — I — l~ l — 1. Applying the conversion by means of 
equations (37) and assuming a spatial discretisation AX = 7.5 m, results 
in a real-world jam density k'- = 1000 4- 7.5 m 133 vehicles/kilometre. 
Conversely, if we know kj, then we can derive kj (see section 2.4.2) and 
hence we have a method to pick a AX. 

If we were to consider multi-cell traffic, e.g., vehicles with different lengths, 
then the jam density would be inversely proportional to the average vehicle 
length. A solution here is to assume a common unit for all vehicle lengths, 
e.g., passenger car units (PCU) [51]. Even though the jam density can be 
defined for each vehicle class separately, it would be more correct to speak of 
an average jam density at this point due to the temporal and spatial variations 
in traffic flows. 



2.4.2 From the real world to a TCA model 

Based on equations (36), we can write the reverse conversion of the microscopic 
characteristics in the following manner: 



h. 



h. 



h' 



9s 



AX 



K 

AT' 



9u = 



AX 



9u 
AT' 



Pi 



AX 



A. 

AT' 



(38) 



In order to agree with our previously stated convention, i.e., all spatial microscopic 
characteristics in a TCA model are integers, equations (38) implicitly assume that 
the real-world spatial variables are multiples of the spatial discretisation (e.g., h' = 
m ■ AX with m G N ). 

Another possible approach to the spatial conversion to TCA model units, is to ap- 
proximate the real-world values as best as possible, whilst keeping our adopted 
convention. As k > 1 cell, this leads to the following scheme where we use up- 
ward rounding (i.e., ceiling): 



\ h '« ] 




' l'i 


AX 


j h 


AX 



(39) 



9 Si h- 
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For example, if AX = 7.5 m, l[ = 4.5 m, and g' s . = 5 m, then h' s . = 4.5+5 = 9.5 m, 
and from equation (39) it follows that h Si = 2 cells, k — I cell, and g Si = 2— 1 = 1 
cell. Because equation (39) is only an approximation, it more than often occurs that 
the computed space headway 'exceeds' the real- world space headway. 

In similar spirit, the conversion for the macroscopic characteristics can be easily 
derived from equations (37). However, as opposed to equations (38) and (39), there 
is no need for an approximation by means of rounding, because these quantities are 
real numbers, as mentioned in the introduction of section 2.4. 



3 Single-cell models 

Having discussed the mathematical and physical aspects of cellular automata and 
TCA models in particular, we now focus on single-cell models. As explained be- 
fore in section 1.3, each cell can either be empty, or is occupied by exactly one 
vehicle; all vehicles have the same length k — I cell. Traffic is also considered to 
be homogeneous, so all vehicles' characteristics are assumed to be the same. In the 
subsequent sections, we take a look at the following TCA models (accompanied by 
their suggested abbreviations): 

• Deterministic models 

■ Wolfram's rule 184 (CA-184) 

• Deterministic Fukui-Ishibashi TCA (DFI-TCA) 

• Stochastic models 

■ Nagel-Schreckenberg TCA (STCA) 

• STCA with cruise control (STCA-CC) 

• Stochastic Fukui-Ishibashi TCA (SFI-TCA) 

■ Totally asymmetric simple exclusion process (TASEP) 

• Emmerich-Rank TCA (ER-TCA) 

• Slow -to -start models 

■ Takayasu-Takayasu TCA (T 2 -TCA) 

■ Benjamin, Johnson, and Hui TCA (BJH-TCA) 

■ Velocity-dependent randomisation TCA (VDR-TCA) 

• Time-oriented TCA (TOCA) 

• TCA models incorporating anticipation 

■ Ultra discretisation, slow-to-accelerate, and driver's perspective 

For other excellent overviews of TCA models, we refer the reader to the works of 
Chowdhury et al. [3], Knospe et al. [19], Nagel [37], Nagel et al. [52], Schadschnei- 
der [53,54], and Schreckenberg et al. [55]. 

All following TCA models will be empirically studied using simulations that are 
performed on a unidirectional, single-lane lattice with periodic boundary condi- 
tions, i.e., a closed loop with Lc — 1. The length of this lattice equals K c = 10 3 
cells, which is taken large enough in order to reduce most unwanted finite-size 
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effects. Our own experiments indicate that larger lattice sizes do not render any sig- 
nificant advantage, aside from the burden of a larger computation time. 




The importance of studying closed-loop, single-lane traffic 

There is often a criticism expressed as to why it is important to study the 
behaviour of traffic flows in such a simplified system. After all, can such a 
basic system capture all the dynamics of real-life traffic flows, or be even 
representative of them ? The answer to this question is that, in our opinion, 
the dynamics of these constrained systems play an important, non-negligible 
role. For example, when considering traffic flows on most unidirectional 
two-lane European motorways, drivers are by law obliged to drive on the 
right shoulder lane, unless when performing overtaking manoeuvres. A fre- 
quently observed phenomenon is then that under light traffic conditions (e.g., 
10 vehicles/km/lane), a slower moving vehicle (e.g., a truck) is located on the 
right lane, and is acting as a moving bottleneck. As a result, all faster vehi- 
cles will line up on the left lane (overtaking on the right lane is prohibited 
by law), thereby causing a density or lane inversion [56,57,5,58]. It is under 
these circumstances that the stability of the car-following behaviour plays an 
important role. Similarly, in densely congested traffic, e.g., the synchronised- 
flow regime, the same stability may govern the fact whether or not a traffic 
breakdown is likely to be induced (see our work in [1] for a discussion on the 
nature of this breakdown). Even for multi-lane traffic, we believe its dynam- 
ics are essentially those of parallel single lanes when considering densely 
congested traffic flows. Another argument for the necessity of studying these 
simplified systems, is the one given by Nagel and Nelson. They state that this 
is the easiest way to determine whether or not internal effects of a traffic flow 
model play a role in e.g., the spontaneous breakdown of traffic, as all exter- 
nal effects (i.e., the boundary conditions) are eliminated [59]. Nevertheless, 
when applying these models to real-life traffic networks, closed-loop traffic 
is not very representative, as the behaviour near bottlenecks plays a far more 
important role [6]. 



All measurements on the TCA models' lattices are based on two possible initial 
conditions: depending on the nature of the study, we will either use homogeneous 
initial conditions (the default), or a compact superjam to start with. In the former 
case, all vehicles are uniformly distributed over the lattice, implying equal space 
headways. In the latter case, all vehicles are 'bunched up' behind each other, with 
zero space gaps. When going from one global density to another, an equivalent 
method would be to adiabatically add (or remove) vehicles to an already homo- 
geneous or jammed state. In our experiments however, we always reset the initial 
conditions, corresponding to the first method. The simulations ran each time for 10 4 
time steps, after an initial period of 10 3 time steps was discarded in order to let tran- 
sients from the initial conditions in the system die out. Global densities, flows, and 
space-mean speeds are computed by means of equations (29) - (32), whereas we 
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use a point detector, i.e., equations (33) - (35), to compute their local variants. In 
this latter case, the data points were collected with a measurement period T mp = 60 
time steps. Based on these results, we can construct (k g ,v Sg ), (k g ,q g ), (ki,v Sl ), and 
(ki,qi) diagrams. To keep a clear formulation, we will however from now on drop 
the subscripts denoting global and local measurements. All simulations were per- 
formed by means of our Traffic Cellular Automata + software (developed for the 
Java™ Virtual Machine); more information can be found in appendix A. 

For a deeper insight into the behaviour of the space-mean speed v s , the average 
space gap g s , and the median time gap g t , detailed histograms showing their distri- 
butions are provided. Note that with respect to the time gaps and time headways, we 
will work in the remainder of this report with the median instead of the arithmetic 
mean. The median gives more robust results when h U) g t . — > +oo, which occurs 
when a vehicle i stops. These histograms are interesting because in the existing 
literature (e.g., [60,53,6]) these distributions are only considered at several distinct 
global densities, whereas we show them for all densities. Each of our histograms is 
constructed by varying the global density k between and 1 , computing the space- 
mean speed, the average space gap and the median time gap for each simulation 
run. A simulation run consists of 5 x 10 4 time steps (with a transient period of 500 
time steps) on systems of 300 cells, varying the density in 150 steps. Note that a 
larger size of the system's lattice, has no significant effects on the results, except 
for an increase of the variance [61]. 




Before giving an elaborate discussion of some of the classic TCA models, it 
is worthwhile to mention the first historical and practical implementations of 
traffic cellular automata. Cremer and Ludwig conceived an implementation 
of traffic flows based on lattice gas automata (LGA), which are a special 
case of cellular automata typically employed when simulating viscous flu- 
ids [62]. Their seminal work, using individual bits to represent vehicles, was 
extended by Schiitt, who provided a simulation package for heterogeneous 
traffic, multi-lane motorways, and network and city traffic [63]. Unfortu- 
nately, the developed models were quite inefficient when they were used in 
setting that called for large scale Monte Carlo simulations [64] . 



3.1 Deterministic models 

In this section, we discuss Wolfram's original rule 184, and its generalisation to 
higher speeds as proposed by Fukui and Ishibashi's deterministic model. We ab- 
breviate these two TCA models as CA- 1 84 and DFI-TCA, respectively. 

3.1.1 Wolfram's rule 184 (CA-184) 

The first deterministic model we consider, is a one-dimensional TCA model with 
binary states. As L c = 1, this model is called an elementary cellular automaton 
(ECA), according to the terminology introduced in section 1.2. If we furthermore 
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assume a local neighbourhood of three cells wide (i.e., a radius of 1), then there are 
2 2 = 256 different rules possible, according to equation (3). Around 1983, Stephen 
Wolfram classified all these 256 binary ECAs [21]. One of these is called rule 184, 
who's name is derived from Wolfram's naming scheme. 

Wolfram's scheme is based on the representation of how a cell's state evolves in 
time, depending on its local neighbourhood. In Fig. 6, we have provided a conve- 
nient visualisation for the evolution of the states in a binary ECA. Here, we can see 
the state Ui{t) of a central cell i at time step t, together with the states (Xi-i(t) and 
(t) of its two direct neighbours i — 1 and i + 1, respectively. All three of them 
constitute the local neighbourhood Mi(t) of radius 1 (see also our example of a CA 
in section 2.1.2). Because states are binary, we can indicate them with a colour, 
i.e., a black square represents a state of 1 (e.g., state <7i+\(t) in Fig. 6), whereas an 
empty (white) square represents a state of 0. According to the local transition rule 
5(i, t), the local neighbourhood A/"j(t) is then mapped from t to t+l onto anew state 
Uiit + 1). The graphical representation in Fig. 6 thus provides us with an illustrative 
method to indicate the evolution of {ovi^), Ui(t), a i+ i(t)} i — > a^t + 1). 

i — li i + 1 
t I I ■ 



t + l 



Fig. 6. An illustrative method for representing the evolution of a cell's state in time, based 
on its local neighbourhood. We can see the state <Ji{t) of a central cell i at time step t, 
together with the states <7j_i(£) and crj+i(i) of its two direct neighbours i — 1 and i + 1, 
respectively. This local neighbourhood is mapped onto a new state <jj(£ + 1). For binary 
states, we use a black square to represent a state of 1 (e.g., state <7j+i(i)), and an empty 
(white) square for a state of 0. The depicted transition maps the triplet (001)2 onto the state 

<7i(t+l) =0. 

Considering the transition depicted in Fig. 6, we can see that a complete neighbour- 
hood contains three cells, each of which can be in a (white) or 1 (black) state. So 
in total, there are 2 3 = 8 possible configurations for such a local neighbourhood. 
Wolfram's naming scheme for the binary ECAs is now based on an integer coding 
of this neighbourhood. Indeed, the local transition rule 5(i,t) is given by a table 
lookup containing eight entries, one for each of the possible local neighbourhoods. 
If we binary sort these eight configurations in the descending order (111), (110), 
(101), (100), (01 1), . . . , then we obtain a graphic scheme such as the one in Fig. 7. 
As can be seen, for each of the local configurations, a resulting or 1 state is re- 
turned for cell % at time step t + l. Collecting all resulting states, and writing them 
in base 2, results in the number (101 1 1000) 2 . Converting this code to base 10, we 
obtain the number 184. Wolfram now coded all 256 possible binary ECAs by a 
unique number in the range from to 255, resulting in 256 rules for these CAs. 

Rule 184 (which we abbreviate as CA-184) is an asymmetrical rule because #((110)2, t) = 
7^ 5((011)2,0 = 1. It is also called a quiescent rule because #((000) 2 ,t) = 
(so all zero-initial conditions remain zero). As an example of the rule's evolution, 
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Ill 110 101 100 011 010 001 000 

10 1110 

Fig. 7. A graphical representation of Wolfram's rule 184, which is written as (10111000) 2 
in base 2. All 8 possible configurations for the local neighbourhood are sorted in descending 
order, expressing the local transition rule S(i,t) as explained by Fig. 6. For example, the 
local neighbourhood (100)2 gets mapped onto a state of 1. This has the physical meaning 
that a particle (black square) moves to the right if its neighbouring cell is empty. 

Fig. 7 shows that the local neighbourhood (100) 2 gets mapped onto a state of 1. 
If we consider these 1 states as particles (i.e., vehicles), and the states as holes, 
then rule 1 84 dictates that all particles move one cell to the right, on the condition 
that this right neighbour cell is empty. Equivalently, all holes have the tendency to 
move to the left for each particle that moves to the right, a phenomenon which is 
termed the particle-hole symmetry. 

For a TCA model, we can rewrite the previous actions as a set of rules that are 
consecutively applied to all vehicles in the lattice, as explained in section 1 .4. For 
the CA-184, we have the following two rules: 

Rl: acceleration and braking 

Vi{t)^wuL{g a .{t-\),\}, (40) 

R2: vehicle movement 

x t (t) ^ x t (t-l)+v t (t). (41) 

Rule Rl, equation (40), sets the speed of the i th vehicle, for the current updated 
configuration of the system; it states that a vehicle always strives to drive at a speed 
of 1 cell/time step, unless its impeded by its direct leader, in which case g Si (t — 1) = 
and the vehicle consequently stops in order to avoid a collision. The second rule 
R2, equation (41), is not actually a 'real' rule; it just allows the vehicles to advance 
in the system. 

In Fig. 8, we have applied these rules to a lattice consisting of 300 cells (closed 
loop), showing the evolution over a period of 580 time steps. The time and space 
axes are oriented from left to right, and bottom to top, respectively. In the left part, 
we show a free-flow regime with a global density k = 0.2 vehicles/cell, in the 
right part we have a congested regime with k = 0.75 vehicles/cell. Each vehicle 
is represented as a single coloured dot; as time advances, vehicles move to the 
upper right corner, whereas congestion waves move to the lower right corner, i.e., 
backwards in space. From both parts of Fig. 8, we can see that the CA-184 TCA 
model constitutes a fully deterministic system that continuously repeats itself. A 
characteristic of the encountered congestion waves is that they have an eternal life 
time in the system. 
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Fig. 8. Typical time-space diagrams of the CA-184 TCA model. The shown closed-loop 
lattices each contain 300 cells, with a visible period of 580 time steps (each vehicle is 
represented as a single coloured dot). Left: vehicles driving a free-flow regime with a global 
density k = 0.2 vehicles/cell. Right: vehicles driving in a congested regime with k = 0.75 
vehicles/cell. The congestion waves can be seen as propagating in the opposite direction of 
traffic; they have an eternal life time in the system. Both time-space diagrams show a fully 
deterministic system that continuously repeats itself. 

In Fig. 9, we have plotted both the (k,v s ) and (k,q) diagrams. As can be seen from 
the left part, the global space-mean speed remains constant at v s = 1 cell/time step, 
until the critical density k c = 0.5 is reached, at which point v s will start to diminish 
towards zero where the critical density kj = 1 is reached. Similarly, the global flow 
first increases and then decreases linearly with the density, below and respectively 
above, the critical density. Here, the capacity flow q cap = 0.5 vehicles/time step is 
reached. The transition from the free-flowing to the congested regime is charac- 
terised by a first-order phase transition. As is evidenced by the isosceles triangular 
shape of the CA-184's resulting (k,q) fundamental diagram, there are only two 
possible kinematic wave speeds, i.e., +1 and -1 cell/time step. Both speeds are also 
clearly visible in the left, respectively right, time- space diagrams of Fig. 8. More 
analytical details on these values will be provided in the following section 3.1.2. 




Global density [vehicles/cell] Global density [vehicles/cell] 



Fig. 9. Left: the (k,v s ) diagram for the CA-184, based on global measurements on the 
lattice. The global space-mean speed remains constant at v s = 1 cell/time step, until the 
critical density k c = 0.5 is reached, at which point v s will start to diminish towards zero. 
Right: the CA-184's (k,q) diagram, with its characteristic isosceles triangular shape. The 
transition between the free-flowing and the congested regime is of a first-order nature. 



26 



3.1.2 Deterministic Fukui-Ishibashi TCA ( DFI- TCA ) 

In 1996, Fukui and Ishibashi constructed a generalisation of the prototypical CA- 
184 TCA model [65]. Although their model is essentially a stochastic one (see sec- 
tion 3.2.3), we will first discuss its deterministic version. Fukui and Ishibashi's idea 
was two-fold: on the one hand, the maximum speed was increased from 1 to v max 
cells/time step, and on the other hand, vehicles would accelerate instantaneously 
to the highest possible speed. Corresponding to the definitions of the rule set of a 
TCA model, the CA-184's rule Rl, equation (40), changes as follows: 

Rl: acceleration and braking 

Vi(t) <- min{g Si (t - (42) 

Just as before, a vehicle will now avoid a collision by taking into account the size of 
its space gap. To this end, it will apply an instantaneous deceleration: for example, 
a fast-moving vehicle might have to come to a complete stop when nearing the end 
of a jam, thereby abruptly dropping its speed from t> max to in one time step. 

Due to the strictly deterministic behaviour of the system, the time-space diagrams 
of the DFI-TCA do not differ much from those of the CA-184. The only difference 
is the speed of the vehicles in the free-flow regime, leading to steeper trajectories. 
It is however interesting to study the (k,v s ) and (k,q) diagrams in Fig. 10. Here 
we can see that increasing the maximum speed v max creates — as expected — a 
steeper free-flow branch in the (k,q) diagram. Interestingly, the slope of the con- 
gested branch does not change, logically implying that the kinematic wave speed 
for jams remains constant, i.e., -1 cell/time step. This can be confirmed with an 
analytical kinematic wave analysis, as explained by Nagel [52]. 




Global density [vehicles/cell] Global density [vehicles/cell] 

Fig. 10. Left: several (k,v s ) diagrams for the deterministic DFI-TCA, each for a different 
^max £ {1, • • • , 5}. Similarly to the CA-184, the global space-mean speed remains constant, 
until the critical density is reached, at which point v s will start to diminish towards zero. 
Right: several of the DFI-TCAs (k,q) diagrams, each having a triangular shape (with the 
slope of the congestion branch invariant for the different v max ). 

Based on the behaviour of the vehicles near the critical density, we can analytically 
compute the capacity flow as follows: in the free-flow regime, all vehicles move 



27 



with a constant speed of v max cells/time step. When the critical density is reached, 
all vehicles drive collision-free at this maximum speed, which implies that g Si = 
^max cells. The space headway h Si = v max + 1 (because k — I for single-cell 
models). Consequently, the value for the critical density as [51]: 



kc = T~ = — "XT' (43) 

ft-Sc ^max ~r t 

The capacity flow is now computed by means of the fundamental relation, i.e., 

Qcap k c V m3x . 



fap = f^. (44) 

u max ~ L 

Applying equations (43) and (44), for e.g., t> max = 5 cells/time step, results in 
k c PS 0.167 vehicles/cell and g cap ps 0.83 vehicles/time step. If we furthermore 
assume AX = 7.5 m and AT = 1 s, then both values correspond to 22 vehi- 
cles/kilometre and 3000 vehicles/hour, respectively. 

As opposed to the instantaneous acceleration in rule Rl, equation (42), we can 
also assume a gradual acceleration of one cell per time step (the braking remains 
instantaneous): 

Rl: acceleration and braking 



Vi(t) <- mm{vi(t - 1) + l,g Si (t - l),v max }. (45) 

However, our experimental observations have indicated that there is no difference 
in global system dynamics, with respect to either adopting gradual or instantaneous 
vehicle accelerations. 




There exists a strong relation between the previously discussed deterministic 
TCA models, and the macroscopic first-order LWR model with a triangular 
q e (k) fundamental diagram [1]. Some of the finer results in this case, are 
the work of Nagel who extensively discusses some analytical results of both 
deterministic and stochastic TCA models [37], and the work of Daganzo who 
explicitly proves an equivalency between two TCA models and the kinematic 
wave model with a triangular q e (k) fundamental diagram [66]. More details 
with respect to such analytical relations, are given in sections 3.2.4 and 5.3. 



To conclude our discussion of deterministic models, we take a look at what hap- 
pens in the limiting case where v max — > +oo. As can be seen in Fig. 11, the con- 
gested branches in both (k,v s ) and (k,q) diagrams grow, at the cost of the free- 
flow branches which disappear. Interestingly, these diagrams correspond one-to- 
one with a triangular q e (k) fundamental diagram that is now expressed in a moving 
coordinate system, as explained by Newell [67]. In such a simplified system, the 
critical density k c = 0, with a capacity flow g cap = 1 . 
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Fig. 11. Left: the (k,v s ) diagram for the deterministic CA-184, with now v max — > +00. 
Right: the (k,q) diagram for the same TCA model, resulting in a critical density k c = 0, 
with a capacity flow q cap = 1 . This type of diagram corresponds to a simplified triangular 
q e (k) fundamental diagram that is expressed in a moving coordinate system. 

3.2 Stochastic models 

The encountered models in the previous section were all deterministic in nature, im- 
plying that there can be no spontaneous formation of jam structures. All congested 
conditions produced in those models, essentially stemmed from the assumed initial 
conditions. In contrast to this, we now discuss stochastic TCA models (i.e., these 
are probabilistic CAs) that allow for the spontaneous emergence of phantom jams. 
As will be shown, all these models explicitly incorporate a stochastic term in their 
equations, in order to accomplish this kind of real-life behaviour [68]. 

3.2.1 Nagel-Schreckenberg TCA (STCA) 

In 1992, Nagel and Schreckenberg proposed a TCA model that was able to repro- 
duce several characteristics of real-life traffic flows, e.g., the spontaneous emer- 
gence of traffic jams [38,64]. Their model is called the NaSch TCA, but is more 
commonly known as the stochastic traffic cellular automaton (STCA). It explic- 
itly includes a stochastic noise term in one of its rules, which we present in the 
same fashion as those of the previously discussed deterministic TCA models. The 
STCA then comprises the following three rules (note that in Nagel and Schrecken- 
berg's original formulation, they decoupled acceleration and braking, resulting in 
four rules): 

Rl: acceleration and braking 

Vi(t) <- min{vi(t -l) + l,g 8 .(t- 1), v max }, (46) 
R2: randomisation 

£(t) <p^ Vi (t) <- max{0,Vi(t) - 1}, (47) 
R3: vehicle movement 
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Xi(t) <— Xi(t - 1) + Vi(t). 



(48) 



Like in both CA-184 and DFI-TCA deterministic TCA models (see sections 3.1.1 
and 3.1 .2), the STCA contains a rule for increasing the speed of a vehicle and brak- 
ing to avoid collisions, i.e., rule Rl, equation (46), as well as rule R3, equation (48), 
for the actual vehicle movement. However, the STCA also contains an additional 
rule R2, equation (47), which introduces stochasticity in the system. At each time 
step t, a random number e [0, 1[ is drawn from a uniform distribution. This 
number is then compared with a stochastic noise parameter p E [0, 1] (called the 
slowdown probability); as a result, there is a probability of p that a vehicle will 
slow down to Vi(t) — 1 cells/time step. The STCA model is called a minimal model, 
in the sense that all these rules are a necessity for mimicking the basic features of 
real-life traffic flows. 

I According to Nagel and Schreckenberg, the randomisation of rule R2 
captures natural speed fluctuations due to human behaviour or vary- 
ing external conditions. The rule introduces overreactions of drivers when 
braking, providing the key to the formation of spontaneously emerging jams. 
Although the above rationale is widely agreed upon, much criticism was 
however expressed due to this second rule. For example, Brilon and Wu be- 
lieve that this rule has no theoretical background and is in fact introduced 
quite heuristically [69]. 

To get an intuitive feeling for the STCA's system dynamics, we have provided two 
time- space diagrams in Fig. 12. Both diagrams show the evolution for a global den- 
sity of k — 0.2 vehicles/cell, but with p set to 0.1 for the left diagram, and p = 0.5 
for the right diagram. As can be seen in both diagrams, the randomisation in the 
model gives rise to many unstable artificial phantom mini-jams. The downstream 
fronts of these jams smear out, forming unstable interfaces [52]. This is a direct 
result of the fact that the intrinsic noise (as embodied by p) in the STCA model 
is too strong: a jam can always form at any density, meaning that breakdown can 
(and will) occur, even in the free-flow traffic regime. For low enough densities how- 
ever, these jams can vanish as they are absorbed by vehicles with sufficient space 
headways, or by new jams in the system [70]. It has been experimentally shown 
that below the critical density, these jams have finite life times with a cut-off that 
is about 5 x 10 5 time steps and independent of the lattice size. When the critical 
density is crossed, these long-lived jams evolve into jams with an infinite life time, 
i.e., they will survive for an infinitely long time [71,64,72]. 

In free-flow traffic, a vehicle's speed will fluctuate between v max and v max — 1 , due to 
the randomisation rule R2. We can compute the space-mean speed in the free-flow 
regime by means of a weighted average. This average corresponds to the probability 
1 — p for driving with the speed v max and the probability p for slowing down to the 
speed v max - 1. As such, we get v Sff = T,WiVi/Hwi = [(1 - p)v max + p(v max - 
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Fig. 12. Typical time-space diagrams of the STCA model (similar setup as for the CA-184 
TCA model in Fig. 8). Both diagrams have a global density of k = 0.2 vehicles/cell. Left: 
the evolution of the system for p =0.1. Right: the evolution of the system, but now for 
p =0.5. The effects of the randomisation rule R2 are clearly visible in both diagrams, as 
there occur many unstable artificial phantom mini-jams. Furthermore, the speed w of the 
backward propagating kinematic waves decreases with an increasing p. 

!)]/[(! ~ P) + p\ = v max — P- In agreement with the space-mean speed observed 
in the left (k,v s ) diagram of Fig. 13, we can state that a vehicle will drive with an 
average free-flow speed of vs = v max — P- 




Global density [vehicles/cell] Global density [vehicles/cell] 



Fig. 13. Left: several (k,v s ) diagrams for the STCA, each for a different p £ {0.1, 0.5, 0.9}. 
It is clear from the diagram, that a vehicle will drive with an average free-flow speed of 
vn = w m ax — V- Right: several (k,q) diagrams for the same STCA models as before. The 
slope of the congested branch tends toward zero for an increasing slowdown probability 
p. Note that the seemingly small capacity drops at the critical density in the right part, are 
actually finite-size effects [73,74]. 

As mentioned in section 3.1.2, the slope of the free-flow branch in a (k,q) diagram 
can be changed by adjusting t> max . Similarly, the slope of the congested branch can 
be changed by tuning the slowdown probability p (note that this also affects the 
average free-flow speed). Looking at the (k,q) diagram in the right part of Fig. 13, 
we note that an increase of p will on the one hand result in a smaller ?J ft , and on the 
other hand the congested branch will lie lower, with a smaller critical density k c . 
In this latter case, the speed w of the backward propagating kinematic waves will 
decrease, an effect that is also visible in the time-space diagrams of Fig. 12. Note 
that the presence of noise in the STCA model causes both free-flow and congested 
branches of the (k,q) diagram to be slightly curved, as opposed to the perfectly 
linear branches of the deterministic models. 



31 



If we setp = 0, then the STCA model becomes deterministic; additionally, setting 
t> max will recover the CA-184 TCA model. In the other deterministic case, when 
p = 1, the system behaves differently: in the congested state, all vehicles will come 
to a full stop, thereby reducing the global flow in the system to zero. As a result, 
the congested branch in the (k,q) regime will coincide with the horizontal axis. This 
implies that the behaviour of a system with v max and p = 1 is totally different than 
that of a system with v max — 1 and p = 0. 




Fig. 14. Three (k,q) diagrams based on local measurements in the STCA model with 
^max = 5 cells/time step. Left: p = 0. Middle: p = I. Right: p = | . Points obtained in the 
free-flow regime (i.e., for v s v max cells/time step) are marked with a o, points obtained 
in the congested regime with a •, and points that imply heavy congestion (i.e., for v s < 1 
cell/time step) with a *. Note that for these local diagrams, the slopes of the congested 
branches (indicated by the points marked as *) are the negative of its corresponding slope 
in a global diagram. 

Considering local measurements of the density, flow, and space-mean speed, the 
(k,q) diagrams in Fig. 14 reveal that an increasing slowdown probability p, results 
in (i) a lower value for the critical density, (ii) a lower capacity flow, and (iii) a more 
localised scatter of the data points. 

In Fig. 15, we have plotted a histogram of the distributions of the STCAs vehicles' 
space gaps, for all global densities k e [0, 1]. For very low densities, the distri- 
butions have a distinct maximum, indicating that all vehicles travel with very large 
space gaps. At higher densities, the maxima of the distributions shift toward smaller 
space gaps, as more and more vehicles encounter jams, even leading to a reduction 
of their space gap to zero. Around the critical density however, the distributions are 
smeared out across consecutive densities, but for each of those densities they ex- 
hibit a bimodal structure. Because the STCA contains many jams, the system now 
contains both vehicles in free-flow traffic, as well as vehicles that are in a congested 
state (i.e., driving closer to each other) [60,4,6,72]. 

In similarly spirit, Fig. 16 shows the distribution of the vehicles' speeds and time 
gaps. Corresponding with our observations of the (k,v s ) diagrams in Fig. 13, the left 
part of Fig. 16 shows a distinct cluster of probability mass at the histogram class 
^max — V f° r ver y iOW global densities. In this region, the standard deviation of the 
space-mean speed is more or less constant and equal to p. At higher global densities, 
the distributions become temporarily bimodal, after which they again tend to a 
unique maximum of cells/time step, corresponding to severely congested traffic; 
the standard deviation drastically encounters a maximum at the critical density, 
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Fig. 15. A histogram of the distributions of the vehicles' space gaps g s , as a function of the 
global density k in the STCA (with v max = 5 cells/time step and p = 0.5). In the contour 
plot to the left, the thick solid line denotes the average space gap, whereas the thin solid 
line shows its standard deviation. The grey regions denote the probability densities. The 
histograms (A) and (B) to the right, show two cross sections made in the left contour plot 
at k = 0.1325 and k = 0.4000 respectively: for example, in (B), the distribution exhibits 
a distinct unique maximum at the histogram class g s = cells, corresponding to the dark 
region in the lower right corner of the contour plot where high global densities occur. 

after which it declines steadily. With respect to the distributions of the time gaps, 
the right part of Fig. 16 shows an rapidly decreasing median time gap as the critical 
density is approached. At this density, the time gaps settle around a local cluster 
at the minimum of 1 time step. Going to higher global densities, the number of 
stopped vehicles increases rapidly, frequently resulting in infinite time gaps. From 
the critical density on, all distributions exhibit a bimodal structure, corresponding 
to vehicles that are caught inside a jam, and other vehicles that are able to move 
freely (possibly at a lower speed) [60,75,72]. 




Fig. 16. Histograms of the distributions of the vehicles' speeds v {left) and time gaps g t 
(right), as a function of the global density k in the STCA (with v max = 5 cells/time step 
and p = 0.5). The thick solid lines denote the space-mean speed and median time gap, 
whereas the thin solid line shows the former's standard deviation. The grey regions denote 
the probability densities. 
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3.2.2 STCA with cruise control (STCA-CC) 

As mentioned in the previous section 3.2.1, a typical artifact of the STCA model is 
that it gives rise to many unstable artificial jams. Due to the noise inherent in the 
model, a jam can always form at any density, even inducing a local breakdown of 
traffic in the free-flow traffic regime. One way to remedy this, is by stabilising the 
free-flow branch of the (k,q) diagram. This can be done by inhibiting the randomi- 
sation for high-speed vehicles. To this end, Nagel and Paczuski considered again 
the rules Rl - R3 of the STCA, i.e., equations (46) - (48), but now complemented 
with a rule RO [73]: 

RO: determine stochastic noise 

Vi(t-l)= v max =>- p'(t) <- 0, ^ 

^ Vi(t - 1) < v max =>> p'(t) <- p, 

with now p replaced by p'(t) in the STCAs randomisation rule R2, i.e., equation 
(47). This new rule effectively turns off the randomisation for high-speed vehicles, 
as only 'jammed' vehicles will now have stochastic behaviour. The resulting TCA 
model, is called the STCA in the cruise-control limit, or STCA-CC for short. If we 
set the maximum speed v max = 1 cell/time step, then all jams initially present in 
the system will coalesce with each other, giving rise to one superjam as depicted in 
Fig. 17. This superjam has been found to follow a random walk in the time-space 
diagram [73,37]. Note that v max > 1 cell/time step does not alter the critical be- 
haviour of the model, even though jam clusters are now branching, having regions 
of free-flow traffic in between them [37]. 




Fig. 17. A time-space diagram of the STCA-CC model for v max = 1 cell/time step and a 
global density of k = 0.4 vehicles/cell. The shown lattice contains 300 cells, with a visible 
period of 1000 time steps. We can see over ten initial jams evolving, coalescing over time 
into one superjam. The system exhibits two distinct phases, i.e., a free-flow and a congested 
regime with v s = 1 and v s = cells/time step respectively. 

In Fig. 18, we show the (k,v s ) and (k,q) diagram of the STCA-CC with v max = 5 
cells/time step and p = 0.2. As can be seen in the right part, the (k,q) diagram 
has a typical inverted A shape (see also our discussion in [1] about the hysteresis 
and capacity drop phenomena). The STCA-CC is said to be bistable, in that both 
the free-flow as well as the congested branches of the (k,q) diagram are stable (the 
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former because it is noise-free). Vehicles going from the free-flow to the congested 
regime encounter at the critical density a phenomenon much like a capacity drop. 
The reverse transition to the free-flow branch proceeds via a lower density and, 
correspondingly, a lower flow (which is the outflow g out of a jam). Comparing the 
right parts of Fig. 13 and Fig. 18, it is evident that a destabilisation of the free-flow 
branch forms the main reason for a lower capacity flow, reached at a lower critical 
density. 




Global density [vehicles/cell] Global density [vehicles/cell] 



Fig. 18. Two (k,v s ) (left) and (k,q) (right) diagrams for the STCA-CC model, with 
^max = 5 cells/time step and p = 0.2. The thick solid line denotes global measurements 
that were obtained when starting from homogeneous initial conditions; the thin solid line 
is based on a compact superjam as the initial condition (see section 3 for an explanation of 
these conditions). The right part clearly shows a typical reversed A shape, which indicates 
a capacity drop. Note that the observed smaller drop in flow for the compact superjam, is 
actually a finite-size effect [73,74]. 

To conclude our discussion of the STCA-CC, we note that the use of cruise control 
as an ADAS can have unintended consequences. The traffic system can be per- 
ceived as having an underlying critical point, at which the life times of jams switch 
from finite to infinite (see our discussion at the beginning of section 3.2.1). The 
existence of this point is closely tied to the self-organised criticality (SOC) of the 
STCA model: the outflow from an infinite jam automatically self-organises to a 
state of maximum attainable flow [76,68,64,77]. Stabilising the free-flow branch 
with cruise-control measures, results on the one hand in traffic higher achievable 
flows which is beneficial, but on the other hand the system is driven closer to its 
critical point which is more dangerous. At this stage, travel times will experience a 
high degree of variability, thereby reducing its predictability [78,73,64]. 

3. 2. 3 Stochastic Fukui-Ishibashi TCA ( SFI- TCA ) 

In section 3.1.2, we discussed the deterministic FI-TCA which is a generalisation 
of the CA-184 TCA model. In their original formulation, Fukui and Ishibashi also 
introduced stochasticity, but now only for vehicles driving at the highest possible 
speed of v max cells/time step [65]. We can express the rules of this model, by con- 
sidering the rules R2 and R3 of the STCA, i.e., equations (47) and (48), but now 
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complemented with the DFI-TCA's rule Rl for instantaneous accelerations, i.e., 
equation (42) of section 3.1.2, and, as in the STCA-CC model, an extra rule RO: 

RO: determine stochastic noise 



with now p replaced by p'(t) in the randomisation rule R2. It can be seen that for 
^max = 1, the SFI-TCA and STCA models are the same. Furthermore, for p = 
the SFI-TCA becomes fully deterministic, and in contrast to the STCA's zero-flow 
behaviour (see section 3.2.1), the SFI-TCAs p — 1 case corresponds to the STCA 
with p = and v max — 1 . 

The rationale behind the specific randomisation in the SFI-TCA model, is that 
drivers who are moving at a high speed, are not able to focus their attention in- 
definitely. As a consequence, there will be fluctuations at these high speeds. As 
such, this corresponds to the opposite of a cruise-control limit, e.g., the STCA-CC 
model. There will be no capacity drop, but the effect on the (k,v s ) diagram is that its 
free-flow branch will become slightly downward curving, starting at v s = v max — p 



To conclude, we mention the related work of Wang et al., who studied the SFI-TCA 
both analytically and numerically, providing an exact result for p = 0, and a close 
approximation for the model with p ^ [79]. Based on the SFI-TCA, Wang et al. 
developed a model that is subtly different. They assumed that drivers do not suffer 
from concentration lapses at high speeds, but are instead only subjected to the ran- 
dom deceleration when they are driving close enough to their direct frontal leaders 
[80]. And finally, we mention the work of Lee et al., who incorporate anticipation 
with respect to a vehicle's changing space gap g s as its leader is driving away. This 
results in a higher capacity flow, as well as the appearance of a synchronised-traffic 
regime, in which vehicles have a lower speed, but are all moving [81]. 

3.2.4 Totally asymmetric simple exclusion process (TASEP) 

The simple exclusion process is a simplified well-known particle transport model 
from non-equilibrium statistical mechanics, defined on a one-dimensional lattice. 
In the case of open boundary conditions (i.e., the bottleneck scenario), particles 
enter the system from the left side at an entry rate a, move through the lattice, 
and leave it at an exit rate (3. The term 'simple exclusion' refers to the fact that a 
cell in the lattice can only be empty, or occupied by one particle. When moving 
through the lattice, particles move one cell to the left with probability 7, and one 
cell to the right with probability 5. When 7 = <5, the process is called the symmetric 
simple exclusion process (SSEP); if 7 7^ 5, then it is called the asymmetric simple 
exclusion process (ASEP) [82]. Finally, if we set 7 = and 5 = 1, the system is 
called the totally asymmetric simple exclusion process (TASEP). If we consider the 
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for k = 0. 
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TASEP as a TCA model, then all vehicles move with t> max = 1 cell/time step to 
their direct right-neighbouring cell, on the condition that this cell is empty. 

Updating the configuration of CA essentially amounts to updating the states of all 
its cells. In general, there are two methods for the update procedure: 

Sequential update 

This updating procedure considers each cell in the lattice one at a time. If all cells 
are considered consecutively, two updating directions are possible: left-to-right and 
right-to-left. There is also a third possibility, called random sequential update. Un- 
der this scheme and with N particles in the lattice, each time step is divided in N 
smaller substeps. At each of these substeps, a random cell (or vehicle) is chosen 
and the CA rules are applied to it. As a consequence of the updating procedure, 
each particle is on average updated after N smaller substeps, which introduces a 
certain amount of noise in the system. We have depicted several typical time-space 
diagrams for the ASEP with 7 = 1 — 5 in Fig. 19. Furthermore note that a hidden 
assumption here is that, after completing a substep, the local information is imme- 
diately available to the whole system, which can violate causality (as information 
is now transmitted through the lattice at an infinite speed). 

Parallel update 

This is the classic update procedure that is used for all TCA models discussed in 
this report. For a parallel update, all cells in the system are updated in one and the 
same time step. Compared to a sequential updating procedure, this one is computa- 
tionally more efficient (note that it is equivalent to a left-to-right sequential update). 
There is however one peculiarity associated with this updating scheme: because all 
particles are considered simultaneously, certain lattice configurations can not exist, 
i.e., the Garden of Eden (GoE) states mentioned in section 1.2. An example of such 
a paradisiacal state, is two vehicles right behind each other, with the following hav- 
ing a non-zero speed. This state would imply that in single-lane traffic, the FIFO 
property was violated and consequently a collision occurred. Such GoE states do 
not exist when using a random sequential update. 

In Fig. 20, we have depicted two time-space diagrams for the TASEP with a random 
sequential updating procedure, operating on a closed loop. As can be seen, the 
diagrams qualitatively look the same, and have some of the same characteristic 
features of the time-space diagrams in Fig. 19. For the TASEP, there is no free-flow 
regime, there are no large jams in the system, and, because of the random sequential 
update, all vehicles continuously have the tendency to collide with each other. As a 
consequence, the system is littered with mini-jams in both the low and high density 
regimes [64,37]. Note that the TASEP with open boundary conditions exhibits a 
very rich behaviour, depending on the values for the entry and exit rates a and (3, 
respectively [83,18,54]. 

With respect to the relations between the TASEP with a random sequential update 
and other models, we mention the following two analogies: on the one hand, the 
LWR first-order macroscopic traffic flow model [1] corresponds to the TASEP in 
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Fig. 19. Typical time-space diagrams of the asymmetric simple exclusion process (ASEP 
model) with a random sequential update and 7 = 1 — 5. The shown lattices each contain 400 
cells, with a visible period of 400 time steps (note that for clarity, the space and time axes 
are located horizontally and vertically, respectively). The global densities in the systems 
were set for each row to k G {0.1,0.3,0.5,0.7,0.9} vehicles/cell. For each column, the 
ASEP's probability to move to the left was set to 7 G {0.1, 0.3, 0.5, 0.7, 0.9}. 




Fig. 20. Typical time-space diagrams of the TASEP model with a random sequential up- 
date. The shown lattices each contains 300 cells, with a visible period of 580 time steps. 
The global density in the system was set to k = 0.3 vehicles/cell (left), and k = 0.7 vehi- 
cles/cell (right). The evolution of the system dynamics qualitatively looks the same in both 
diagrams: the system is littered with mini-jams in both the low and high density regimes. 
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the hydrodynamic limit to a noisy and diffusive conservation law, which can be 
reduced to the LWR model [64,37]. On the other hand, the TASEP corresponds to 
the STCA (see section 3.2.1), but now with v max = 1 cell/time step [3,6]. 

To gain more insight into the macroscopic behaviour of the TASEP with random 
sequential update, we provide its (k,v s ) and (k,q) diagrams in Fig. 21. Looking 
at the (k,v s ) diagram on the left part, we notice that the TASEP with t> max = 1 
cell/time step corresponds exactly to Greenshields' original linear relation between 
the density and the mean speed [84,51]. This in fact is a further testimony of the 
close link between the TASEP and the LWR model with a triangular q e (k) funda- 
mental diagram. Increasing the TASEP's maximum speed, leads to a more curved 
relation, intersection the vertical axis at the point (0,t> max ). In any case, the (k,v s ) 
diagram also reveals the absence of a distinct free-flow branch, corresponding to 
the observations of the large amount of mini-jams for all global densities, as could 
be seen in the time-space diagrams of Fig. 20. 

Studying the (k,q) diagram in the right part of Fig. 21, we can see that the TASEP 
corresponds with the STCA for v max = 1 and an arbitrary slowdown probability 
(e.g., p = 0.1). The diagram also shows how the CA-184 leads to a sharp transi- 
tion between the free-flow and the congested regime, as opposed to the rounded 
peak of capacity flow at k = 0.5 vehicles/cell for the STCA. However, whereas 
the TASEP also has its capacity flow at the same value, there does not occur such 
a phase transition as in the other models. Finally, we can see that increasing the 
maximum speed v max for the TASEP introduces no significant qualitative changes, 
except for a skewing towards lower densities [64]. 

I Note that with respect to the computational complexity of the implemented 
TCA models, most measurements in this report took a few hours to obtain, 
using an Intel P4 2.8 GHz with 512 MB RAM, running the Java™ JDK 1.3.1 
under Windows 2000. In sharp contrast to this, are the computations for the 
TASEP model, which took nearly two weeks to complete. 

3.2.5 Emmerich-Rank TCA (ER-TCA ) 

Whereas the classical STCA model provided a reasonable qualitative agreement 
with real-world observations, Emmerich and Rank addressed the quantitative dis- 
crepancies between the model and real-world data. To this end, they proposed a 
variation on the STCA, extending the influence of the space gap on a vehicles up- 
dated speed [85]. 

In their work, Emmerich and Rank fundamentally modified the STCA in two steps: 
(i) they changed the parallel update procedure to a right-to-left sequential update 
procedure (see section 3.2.4 for more details), and (ii) they changed the behaviour 
of vehicles that are slowing down. In a nutshell, (i) leads to the important result that 
vehicles are now able to drive directly behind each other (i.e., with a zero space gap) 
at high speeds, because the gaps in a traffic stream are used more efficiently. The 
reason is that due to the specific sequential update, a downstream vehicle is moved 
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Fig. 21. A comparison of the (k,v s ) (left) and (k,q) (right) diagrams for the CA-184 with 
^max = 1 (A), the STCA with v max = 1 and p = 0. 1 (*), the TASEP with random sequential 
update and v max = 5 (o), and the TASEP with random sequential update and v maK = 1 
(thick solid line). Some distinct characteristics of the TASEP are the absence of a free-flow 
regime, and for v max = 1 cells/time step, the exact correspondence with Greenshields linear 
relation between the density and the mean speed. 

first (for a closed loop, the vehicle with the largest space gap is chosen first), after 
which the next vehicle upstream will see a larger space gap. 

Just as the STCA can be seen as a special case of the optimal velocity model 
(OVM), based on a linear optimal velocity function (for a description of the OVM, 
we refer the reader to our overview in [1]), the ER-TCA model generalises this 
function by making a vehicle's speed dependent on a variable safe distance and its 
current speed [3]. This affects (ii), i.e., vehicles that are slowing down: when deter- 
mining the new speed of a vehicle, the ER-TCA model first checks if the vehicle is 
within 10 cells of its direct frontal leader. If this is the case, then the vehicle will 
slow down according to a table lookup in a gap-speed matrix M 9s . m . This matrix 
is constructed in such a way that collisions are avoided (i.e., M i:j < min{i, j}): 



M 



(51) 



00000000000 
1111111111 

01222222222 
01233333333 
01234444444 
^0 1234444445/ 

The matrix in equation (51), conveys the idea that lower speeds require lower space 
gaps, and that vehicles tend to keep larger space gaps when travelling at higher 
speeds. This latter effect is also visible in the distribution of the vehicles' space 
gaps, as visualised in the histograms in the left part of Fig. 22, where, in contrast 
to the STCA's space gaps distribution of Fig. 15, large space gaps are observed for 
densities near the critical density. Furthermore, because of this mechanism, vehicles 
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will have smoother decelerations, instead of the abrupt slowing down in the STCA 
model and some of its variations. 




- 



Fig. 22. Histograms of the distributions of the vehicles' space gaps g s (left) and time gaps 
gt (right), as a function of the global density k in the ER-STCA (with v max = 5 cells/time 
step and p = 0.35). The thick solid lines denote the mean space gap and median time gap, 
whereas the thin solid line shows the former's standard deviation. The grey regions denote 
the probability densities. 

To understand some of the system dynamics of the ER-TCA model, we have pro- 
vided several (k,v s ) and (k,q) diagrams in Fig. 23. For p = 0.35, we can see in the 
(k,q) diagram in right part, that the free-flow branch gets curved, implying that ve- 
hicles travel at a slightly lower speed when they approach the capacity-flow regime. 
Because vehicles can travel at high speeds in dense platoons, the ER-TCA model 
can achieve very high capacity flows, even leading to q > 1 vehicle/time step. In 
order to constrain these flows to realistic values, the ER-TCA model needs a quite 
high slowdown probability, e.g., p = 0.35. 

These two effects, i.e., a curving of the free-flow branch and an increased capacity 
flow, are basically what the ER-TCA model is all about, as there is no qualitative 
change in the congested branch of the (k,q) diagram. There are however some se- 
rious drawbacks to the ER-TCA model. First and foremost, the (k,q) diagram is no 
longer non-monotonic for low densities when the sequential update is replaced by 
a parallel one [3,19]. Secondly, the model exhibits too large time headways in the 
free-flow regime when compared with real-world data. This effect is also visible 
in the distribution of the vehicles' time gaps, as depicted in the histograms in the 
right part of Fig. 22, where, in contrast to the STCA's time gaps distribution of 
Fig. 16, a large amount of finite time gaps extends well into the region of medium 
densities. Third, due to the sequential update, the ER-TCA model's downstream 
jam dynamics are unstable, just as in the STCA model [19]. Fourth, as can be seen 
from the (k,v s ) diagram in the left part of Fig. 23, for small slowdown probabilities 
p, the resulting space-mean speed in the system is very unrealistic, even including 
plateaus of constant speed in the congested regime, e.g., the curve associated with 
p = 0.1 (we consider p = as a degenerate case). 
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Fig. 23. Left: several (k,v s ) diagrams for the ER-TCA, each for a different slowdown 
probability p. It is clear from the diagram, that for low values of p, the resulting diagrams 
are unrealistic, including plateaus of constant space-mean speed in the congested regime. 
Right: several (k,q) diagrams for the same ER-TCA models as before. Due to the system 
dynamics in the ER-TCA, very high capacity flows are possible. To constrain these flows, 
the slowdown probability p has to be quite large in order to obtain realistic results. In both 
parts of the figure, the thick solid line denotes the original model of Emmerich and Rank, 
who used a value p = 0.35 as their best fit to experimental data. 

3.3 Slow-to-start models 

In order to obtain a correct behavioural picture of traffic flow breakdown and stable 
jam, it is necessary that a vehicle's minimum time headway or reaction time should 
be smaller than its escape time from a jam, or equivalently, the outflow from a jam 
(i.e., the queue discharge rate) must be lower than its inflow [86,70,87,88,89,52]. 
If this is not the case, as in e.g., the STCA model where both times are exactly 
the same, then all jams will be unstable, as can be seen in the time-space diagram 
of Fig. 12. Because of their unstable jamming behaviour, the previously discussed 
stochastic models, experience neither a capacity drop nor a hysteresis loop, for 
which stable jams are a necessary prerequisite. Although the STCA-CC seems to 
be an exception to this rule, the downstream fronts of its jams are still too unstable, 
in the sense that new jams can emerge all too easily, which is unrealistic behaviour 
with respect to real-life traffic flows [5]. 

As just mentioned, one mechanism that deals with this, is by leaving free-flow 
traffic undisturbed, and by significantly reducing the outflow from a jam once a 
breakdown occurs, thereby stabilising the downstream front of a jam. Instead of 
just eliminating the noise in free-flow traffic in the STCA-CC, this reduced outflow 
can also be accomplished more intuitively, by making the vehicles wait a short 
while longer before accelerating again from stand still. As such, they are said to be 
"slow to start" . 

Note that there exists yet another mechanism that allows for the reproduction of the 
capacity drop and hysteresis phenomena (we will only briefly mention it here). The 
approach followed by Werth, is based on the premise that drivers take into account 
the speed difference with their direct frontal leader, instead of just the space gap 
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as was previously assumed. This leads to Galilei invariant vehicle- vehicle interac- 
tions (i.e., the system dynamics remain the same if a new linear moving coordinate 
system is substituted in the equations). Interestingly, the metastability in this model 
is not due to cruise control or slow-to-start rules, but rather a result of the anticipa- 
tion adopted. The model can exhibit stable dense platoons of fast vehicles, resulting 
in a stabilisation of the free-flow branch, and consequently leading to hysteretic be- 
haviour [90,5,3]. 

With respect to real-world units, we give some typical values associated 
with the capacity drop and hysteresis phenomena (based on [5]): an out- 
flow g out PS 1800 vehicles/hour/lane at an associated density of fc out ps 20 
vehicles/km/lane, with q cap , fc crit , and equal to 2700 vehicles/hour/lane, 
20 vehicles/km/lane, and 140 vehicles/km/lane, respectively. 

3.3.1 Takayasu-Takayasu TCA (T 2 -TCA) 

In 1993, Takayasu and Takayasu proposed a deterministic TCA model, based on the 
CA-184 (see section 3.1.1), that incorporated a delay in acceleration for stopped 
vehicles [91]. Their motivation stems from the fact that high-speed vehicles are in 
general able to decelerate very quickly, but conversely, it takes them a lot longer to 
attain this high speed when they start from a stopped condition. As such, Takayasu 
and Takayasu introduced a delay, based on the rationale that a vehicle will only 
start to move when it recognises movement of its direct frontal leader. Translating 
this into a rule set, we can write the T 2 -TCAs rules based on those of the CA-184, 
but now with the following modifications (note that v max = 1 cell/time step): 

Rl: braking 

Vi(t-1) >g Si (t-l) =}► Vi(t)<-g ai (t-l), (52) 
R2: delayed acceleration 

Vt (t-l) = A <7 ai (i-l)>2 
Viit) <- 1, 

R3: vehicle movement 



Xi(t) <-Xi(t-l)+Vi(t). (54) 

From this rule set it follows that a vehicle will always drive at a speed of one 
cell/time step, unless it has to brake and stop according to rule Rl, equation 52. 
Furthermore, the vehicle is only allowed to accelerate again to this speed of one 
cell/time step, on the condition that it has a sufficiently large space gap in front, 
as dictated by rule R2, equation 53. As a result, the introduced delay is spatial in 
nature, and it only affects stopped vehicles. 

In Fig. 24, we have depicted the resulting (k,v s ) and (k,q) diagrams for the T 2 -TCA 
model. The observed behaviour is similar to that of the STCA-CC model in section 
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3.2.2, in that the T 2 -TCA model also exhibits bistability. Starting from homoge- 
neous initial conditions, the space-mean speed in the system undergoes a sharp 
drop once a vehicle has to stop. The reverse process, i.e., going from the congested 
to free-flow regime, is accompanied by a smooth continuous transition. Takayasu 
and Takayasu state that this corresponds to a second-order phase transition, be- 
cause their order parameter (the sum of the jamming times) follows a power-law 
distribution, with jam times tending to infinity once the system goes beyond the 
critical density. With respect to the T 2 -TCA's tempo-spatial behaviour, we note that 
the critical density for the former transition is located at k c = 0.5 vehicles/cell, at 
which point all vehicles travel at a speed of one cell/time step with all space gaps 
equal to one cell. The density at which the recovery associated with latter transition 
occurs, is equal to k — | vehicles/cell, at which point all vehicles travel at a speed 
of one cell/time step, but now with all space gaps equal to two cells. Fukui and 
Ishibashi later modified the delaying process, resulting in a system that always re- 
laxes to a state in which the space-mean speed oscillates between two values, both 
smaller than one cell/time step [92]. 
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Fig. 24. Two (k,v s ) (left) and (k,q) (right) diagrams for the T 2 -TCA model, with -u max = 1 
cells/time step. The thick solid line denotes global measurements that were obtained when 
starting from homogeneous initial conditions; the thin solid line is based on a compact 
superjam as the initial condition (see section 3 for an explanation of these conditions). The 
right part clearly shows a typical reversed A shape, which indicates a capacity drop. 




The original background for Takayasu and Takayasu's work, was based on 
the presence of so-called II f noise (also known as pink noise or flicker 
noise) in the Fourier transformed density fluctuations of motorway traffic. 
The seemingly random stop-and-go motions of jammed vehicles, could indi- 
cate a chaotic behaviour (as opposed to just statistical noise), closely coupled 
with self-organised criticality (see also the end of section 3.2.2) [68]. In the 
free-flow regime of the T 2 -TCA model, jams have a finite life time leading 
to a flat spectrum, as opposed to the congested regime where jams have an 
infinite life time, leading to a II f spectrum [91]. 



Schadschneider and Schreckenberg later provided a generalisation of the T 2 -TCA 
model: keeping v max = 1 cell/time step, they now modified the braking and acceler- 
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ation behaviour of a vehicle. On the one hand, they kept Takayasu and Takayasu's 
original acceleration rule R2, equation (53), and on the other hand, they allowed 
a vehicle with a space gap of just one cell to accelerate with a slow-to-start prob- 
ability 1 — p t [93]. They furthermore also introduced a randomisation for moving 
vehicles, similar to the STCA (see section 3.2.1), making vehicles stop with a slow- 
down probability p. Several interesting phenomena occur for certain values of both 
probabilities p and p t . The modified spatial slow-to-start rule can lead to the appear- 
ance of an inflection point in the (k,q) diagram at very high densities. The effect gets 
strongly exaggerated when p t — > 1, at which point a completely blocked state of 
zero flow appears for all global densities k > 0.5 vehicles/cell [93,18,3]. 



3.3.2 The model of Benjamin, Johnson, and Hui (BJH-TCA) 

Around the same time that Takayasu and Takayasu proposed their T 2 -TCA model, 
Benjamin, Johnson, and Hui (BJH) constructed another type of TCA model, using 
a slow-to-start rule that is of a temporal nature [94]. Their BJH-TCA model is 
based on the STCA (see section 3.2.1), but extended it with a rule that adds a 
small delays to a stopped car that is pulling away from the downstream front of a 
queue. Benjamin et al. attribute this rule to the fact that it mimics the behaviour 
of a driver who momentarily looses attention, or when a vehicle's engine is slow 
to react. Their slow-to-start rule allows a stopped vehicle to move again with this 
slow-to-start probability 1 — p s . If the vehicle did not move, then it tries to move 
again but this time with probability p s . Due to this peculiar acceleration procedure, 
all vehicles require a memory that, as mentioned before, makes the slow-start-rule 
temporal in nature [3]. As a result of this new systematic behaviour, jams will now 
become less ravelled (as opposed to the STCA), because the slow-to-start rule will 
have the tendency to merge queues. 

The BJH-TCA model was also applied to the description of a motorway with an 
on-ramp, leading to the conclusions that (i) it actually is beneficial to have jams on 
the main motorway, due to the fact that these jams homogenise the traffic streams 
as they compete for stopped vehicles, and (ii) it is desirable to set a maximum 
speed limit on this main motorway which allows to maximise the performance of 
the on-ramp. Note that in their discussion, Benjamin et al. used the queue length 
at the on-ramp as a performance measure. In our opinion, this is not a very good 
choice as it ignores e.g., the total time spent in the system, which we believe is a 
more important measure (see also the work of Bellemans [95] and Hegyi [96] in 
this respect). 

To conclude, we note that the (k,q) diagrams of the BJH-TCA and T 2 -TCA models 
qualitatively look the same, with the exception the former does not have the possi- 
bility of an inflection point, or a density region with zero flow, as was the case for 
the latter model (see section 3.3.1) [93,18]. 
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3.3.3 Velocity-dependent randomisation TCA (VDR-TCA) 

As already explained in the introduction of this section, reducing the outflow from 
a jam is responsible for the capacity drop and hysteresis phenomenon. To this end, 
Barlovic et al. proposed a TCA model that generalises the STCA model (see section 
3.2.1) by employing an intuitive slow-to-start rule for stopped vehicles [97,98]. 
Similar to the STCA-CC (see section 3.2.2), the complete rule set for the VDR- 
TCA is as follows: 

RO: determine stochastic noise 

Vi(t- 1) =0 
Vi{t-\) >0 
Rl: acceleration and braking 

Vi(t) <- mm{vi{t -l) + l,g s .(t- l),v max }, (56) 
R2: randomisation 

£(t) < p'(t) =s> Vi (t) <- max{0, Vi (t) - 1}, (57) 
R3: vehicle movement 

Xi(t)<-Xi(t-1)+Vi(t). (58) 

As before, in rule R2, equation (57), £(£) e [0, 1 [ denotes a uniform random num- 
ber (specifically drawn for vehicle i at time t) and p'(t) is the stochastic noise 
parameter, dependent on the vehicle's speed (hence the name 'velocity-dependent 
randomisation'). The probabilities p and p are called the slow-to-start probability 
and the slowdown probability, respectively, with p ,p e [0, 1]. Note that Barlovic 
et al. only considered the case with two different noise parameters (i.e., p Q and p), 
ignoring the more general case where we can have a noise parameter for each pos- 
sible speed (i.e., p , p Vms J. Their model was also considered for systems with 
open boundary conditions [99] . 

Depending on their speed, vehicles are subject to different randomisations: typical 
metastable behaviour results when p 3> p, meaning that stopped vehicles have to 
wait longer before they can continue their journey. This has the effect of a reduced 
outflow from a jam, so that, in a closed system, this leads to an equilibrium and 
the formation of a compact jam. For such a typical situation, e.g., p = 0.5 and 
p = 0.01, the tempo-spatial evolution is depicted in Fig. 25. We can see an initially 
homogeneous traffic pattern (one metastable phase) breaking down and kicking the 
system into a phase-separated state, consisting of a compact jam surrounded by 
free-flow traffic. In such a state, traffic jams in the system will absorb as many 
vehicles as is necessary, in order to have a free-flow phase in the rest of the system 
[6]. Note that the VDR-TCA can also be equipped with a cruise control, by turning 
of fluctuations for vehicles driving at the maximum speed t> max . 



P'(t) 
p'(t) 



Po, 
P, 



(55) 
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Fig. 25. A time-space diagram of the VDR-TCA model for v max = 5 cells/time step, 
po = 0.5, p = 0.01, and a global density of k = 4 vehicles/cell. The shown lattice contains 
300 cells, with a visible period of 1000 time steps. We can see the breakdown of an initially 
homogeneous traffic pattern. As the phase separation takes place, a persistent compact jam 
is formed, surrounded by free-flow traffic. The significant decrease of the density in the 
regions outside the jam results from the jam's reduced outflow. 

In the left part of Fig. 26, we have plotted a histogram of the distributions of the ve- 
hicles' speeds, for all global densities k e [0, 1]. Here we can clearly see the distinc- 
tion between the free-flow and the congested regime: the space-mean speed remains 
more or less constant at a high value, then encounters a sharp transition (i.e., the ca- 
pacity drop), resulting in a steady declination as the global density increases. Note 
that as the critical density is encountered, the standard deviation jumps steeply; this 
means that vehicles' speeds fluctuate wildly at the transition point (because they 
are entering and exiting the congestion waves). Once the compact jam is formed, 
the dominating speed quickly becomes zero (because vehicles are standing still in- 
side the jam). Although most of the weight is attributed to this zero-speed, there is 
a non-negligible maximum speed present for intermediate densities. If the global 
density is increased further towards the jam density, this maximum speed disap- 
pears and the system settles into a state in which all vehicles either have speed zero 
or one (i.e., systemwide stop-and-go traffic). 

Studying the (k,q) diagram in the right part of Fig. 26, gives us another view of 
this phase transition. We can see a capacity drop taking place at the critical density, 
where traffic in its vicinity behaves in a metastable manner. This metastability is 
characterised by the fact that sufficiently large disturbances of the fragile equilib- 
rium can cause the flow to undergo a sudden decrease, corresponding to a first-order 
phase transition. The state of very high flow is then destroyed and the system settles 
into a phase separated state with a large megajam and a free-flow zone. The large 
jam will persist as long as the density is not significantly lowered, thus implying 
that recovery of traffic from congestion follows a hysteresis loop. In contrast to the 
STCA-CC's bistability, the VDR-TCA model is truly metastable, because now the 
free-flow branch in the (k,q) diagram becomes unstable for large enough perturba- 
tions. Furthermore, the spontaneous formation of jams in the downstream front that 
troubled the STCA, is suppressed in the VDR-TCA model. 
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Fig. 26. Left: a contour plot containing the histograms of the distributions of the vehicles' 
speeds v as a function of the global density k in the VDR-TCA (with v max = 5 cells/time 
step, po = 0.5 and p = 0.01). The thick solid line denotes the space-mean speed, whereas 
the thin solid line shows its standard deviation. The grey regions denote the probability 
densities. Right: a (k,q) diagrams for the same TCA model. The dotted line denotes global 
measurements that were obtained when starting from homogeneous initial conditions; the 
solid line is based on a compact superjam as the initial condition. The right part clearly 
shows a typical reversed A shape, which indicates a capacity drop. 




Note that if p <C p, then the behaviour of the system will be drastically dif- 
ferent. Four distinct traffic regimes emerge in the limiting case where p = 
and p — 1; in this case, the model is called fast-to-start [100]. In these four 
regimes, moving vehicles can never increase their speed once the system has 
settled into an equilibrium. Furthermore, there exists a regime which expe- 
riences forward propagating density waves, corresponding to a non-concave 
region in the system's flow-density relation. For more information, we refer 
to our work in [101] and [102]. 



3. 3. 4 Time-oriented TCA (TOCA) 

Considering the STCA model (see section 3.2.1), Brilon and Wu acknowledged the 
fact that it is quite capable of reproducing traffic dynamics in urban street networks. 
However, they also recognised the fact that the model performed rather inadequate 
when it comes to correctly describing the characteristics of traffic flows on motor- 
ways, e.g., compared to field data of a German motorway. Brilon and Wu blamed 
the unrealistic car-following behaviour of the STCA model for its inferior capabil- 
ities. At the core of their argument, they attributed this to the fact that the STCA 
model is exclusively based on spatial variables (e.g., space headways). In order to 
alleviate these problems, they proposed to use a model that was based on temporal 
variables (e.g., time headways), leading to more realistic vehicle- vehicle interac- 
tions [69]. The rule set for this time-oriented TCA model (TOCA) is as follows: 

Rl: acceleration 
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g Si (t-l)>(v i (t-l)-g ts ) A 

< Pace 

=^ Vi(t) <- min{^(£ - 1) + \,v max }, 



(59) 



R2: braking 



Vi(t) <- min{^(t),^(t - 1)}, 



(60) 



R3: randomisation 



g ai (t-l)<(vi(t-l)-g ta ) A 

&(*) < Pdec 

=>• <- max{^(t) -1,0}, 



(61) 



R4: vehicle movement 



Xi(t) 



Xi(t- l)+Vi(t). 



(62) 



In the above rules, G [0, 1 [ are random numbers drawn from a uniform 

distribution, ~g ts >A.T * s me sa / e time gap, p acc is the acceleration probability, and 
Pdec is the deceleration probability. Because all interactions between vehicles in 
the STCA are bounded by the update time step, their speeds will never oscillate, 
leading to a rigid and stable system. As a consequence of the TOCA's temporal 
rules however, vehicles will now behave more elastically, taking a safe time gap 
into account that allows them to adapt their speeds with a relaxation. In this case, a 
vehicle will resort to emergency braking (i.e., an instantaneous deceleration) only 
if it gets too close to its direct frontal leader [52]. Typical parameter values for the 
TOCA are g t =1.2 time steps and p acc = p dec = 0.9. Brilon and Wu also extended 
their model with rudimentary rules that allowed for lane changes on unidirectional 
multi-lane roads. 

In the left part of Fig. 27, we can see a similar tempo- spatial behaviour as with 
the VDR-TCA (see section 3.3.3), in that an initially homogeneous traffic pattern 
breaks down, resulting in dilute jam that is surrounded by free-flow traffic. The 
major difference between jamming in the VDR-TCA and TOCA models however, 
is that in the former model, vehicles come to a complete stop when entering a jam 
(see Fig. 25). They remain stationary until they can leave the downstream front of 
the queue. In contrast to this, the jams in the TOCA model contain moving vehicles. 
Pushing the global density even further to k — 0.5 vehicles/cell as was done in the 
right part of Fig. 27, results in a fully developed jam that dominates the entire 
system and contains temporarily stopped vehicles. 

Fig. 28 depicts two groups of (k,q) diagrams for the TOCA model, with v max = 5 
cells/time step. The left part shows four diagrams for different combinations of p acc 



andp dec G {(0.9, 0.1), (0.9, 0.9), (0.1, 0.1), (0.1, 0.9)}, each time with g ts =1.2 
time steps. As can be seen, the default case with p acc = p dec = 0.9 leads to an 
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Fig. 27. Typical time-space diagrams of the TOCA model for -u max = 5 cells/time step, 
g ta =1.2 time steps, and p acc = p<jec = 0-9- The global density was set to k = I vehi- 
cles/cell (left) and k = 0.5 vehicles/cell (right). The shown lattices each contain 300 cells, 
with a visible period of 580 time steps. In the left part, we can see the breakdown of an 
initially homogeneous traffic pattern, resulting in dilute jam that is surrounded by free-flow 
traffic. In the right part, we see a fully developed jam, dominating the entire system. As 
can be seen, for moderately light densities, the jams in the TOCA model contain moving 
vehicles. 



inflection point at a moderately high density of k = 0.5 vehicles/cell, resulting in 
two different slopes for the congested branch of the TOCA's (k,q) diagram. At this 
point, vehicles will have average space gaps less than one cell, and because p& ec is 
rather high, vehicles will have the tendency to slow down (and p acc is smaller then 
one, so their acceleration is somewhat inhibited). As a result, a large jam, com- 
parable to the system's size, will dominate tempo-spatial evolution. Furthermore, 
the acceleration probability p acc should take on rather high values, otherwise the 
global flow in the system is too low because vehicles are not accelerating anymore. 
In the right part of Fig. 28, we have shown a large amount of diagrams for differ- 
ent g ts with p acc = p dec = 0.9. Here we can see that, for ~g ts < AT, the resulting 
density-flow curves are non-mono tonic. Higher values for ~g ts in more vehicles that 
drive more cautiously, apparently leading to higher values for the critical density 
and the capacity flow. Note that the seemingly small capacity drops at the end of 
each free-flow branch are in fact finite-size effects [73,74]. 

In their original paper, Brilon and Wu claim that their TOCA model results in a bet- 
ter agreement with empirical data, a fact which is based on a qualitative comparison 
of the (q,v s ) diagrams [69]. Note that, after personal communication with the au- 
thors, it seems they performed a minimisation of the square errors in the (k,v s ) 
diagram. However, in order to get the correct values for calibrating the TOCA's 
parameters, they just manually guessed, without performing a thorough numeri- 
cal optimisation. Despite this optimistic view, Knospe et al. later investigated the 
TOCA model's capabilities more thoroughly. Their conclusions state that a quite 
large value for the deceleration probability p dec is necessary in order to obtain real- 
istic capacity flows. Although the time headway distribution of a jam's downstream 
front in the TOCA model is correct with respect to real-life observations, its down- 
stream front moves too fast due to the large deceleration probability. As a result, 
the jams in the TOCA model are more dilute, as could be seen in Fig. 27 [19]. 
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Fig. 28. Two groups of (k,q) diagrams for the TOCA model, with v max = 5 cells/time step. 
Left: four diagrams for different combinations of p acc and p^, with g ta =1.2 time steps. 
Right: a large amount of diagrams for different g ts with p acc = p^ ec = 0.9. For g ts < AT, 
the resulting density-flow curves are non-monotonic. Note that the seemingly small capac- 
ity drops at the end of each free-flow branch are in fact finite-size effects [73,74]. 

3.3.5 TCA models incorporating anticipation 

One of the models related to anticipative driving (i.e., only taking a leaders' reac- 
tions into account, without predicting them), can be found in the work of KrauB et 
al., who derived a collision-free model based on the STCA (see section 3.2.1), but 
which uses continuous vehicle speeds. Their model can be considered as a sim- 
plified version of the Gipps model [1]. Although the model restricts vehicles' de- 
celeration capabilities, it is still able to correctly reproduce the capacity drop and 
hysteresis phenomena [74]. 

Another model with anticipation was proposed by Eissfeldt and Wagner [41]. Their 
model is based on KrauB's work [1], and employs a next-nearest-neighbour interac- 
tion, which stabilises dense flows and results in a non-unique flow-density relation. 

Recently, Larraga et al. introduced a TCA model that includes a driver's anticipa- 
tion of the leading vehicle's speed [42]. In contrast to the STCA model (see section 
3.2. 1), the acceleration and braking rules are decoupled. As a first rule, the standard 
acceleration towards the maximum speed is applied, after which the randomisation 
is performed by means of a second rule. Only then, the model considers braking in 
its third rule; however, the deceleration is not only based on the space gap between 
both vehicles, but also on an anticipation of the leading vehicle's speed: 

R3: anticipation and braking 



Vi(t) 



mm < 



Vi(t),g Si (t-l) + 



1 - (Xi) ■ v i+ i(t - 1) + - 



(63) 



safe distance 

with Vi (t) on the right-hand side corresponding to the computed speed after apply- 
ing rule R2, [x] denoting x rounded to the nearest integer, v i+ \(t — 1) the speed of 
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the leading vehicle at the current time step, and ccj £ [0, 1] an anticipatory driving 
parameter for the i th vehicle. In their work, Larraga et al. considered all a { to be 
equal. 

The interesting aspect of this anticipatory TCA model, is that for certain values 
of a, it can result in dense platoons of vehicles, travelling coherently and thereby 
leading to forward propagating density structures. In the free-flow regime, the (k,q) 
diagram also exhibits a slight curvature near the capacity flow, similar to the ER- 
TCA model (see section 3.2.5). Del Rio and Larraga later also extended the model 
to accommodate for multi-lane traffic flows [103]. 

3.3.6 Ultra discretisation, slow-to-accelerate, and driver's perspective 

It is also possible to derive a cellular automaton model, based on the discretisation 
of a partial differential equation. Starting from a PDE (e.g., the Burgers equation 
[1], we can obtain an finite difference equation by discretising the spatial and tem- 
poral dimensions, resulting in a model that still has continuous state variables. As 
a further step, we can now also discretise these state variables, using a process 
called the ultra-discretisation method (UDM) [104]. The result of the UDM can 
be interpreted as a cellular automaton in the Euler representation. The latter means 
that for a TCA model, a road is considered to be a field, whereby the individual 
cars are not distinguished [105]. The interesting part of this type of CA is that its 
cells are allowed to hold multiple vehicles, which makes it possible to implicitly 
model multi-lane traffic in a simplified sense (because the effects of lane changes 
are neglected) [3]. As a next step, this obtained CA can be cast in its Lagrangian 
representation, by means of an Euler-Lagrange transformation [105,106]. The re- 
sulting Lagrange representation treats the positions of all vehicles individually, thus 
leading to the well-know position-based rule sets of the TCA models discussed in 
this report. 

Nishinari proposed an interesting TCA model, based on the above UDM scheme. 
Their discretisation leads to the so-called Burgers cellular automaton (BCA), which 
is for single-lane traffic equivalent to the CA-184 TCA model (see section 3.1.1) 
[107,108]. Emmerich et al. also provided a TCA model, by applying the UDM 
scheme to a Korteweg-de Vries equation. In contrast to the BCA model, their work 
resulted in a second-order TCA model because the CA's global map not only needs 
the configuration at the previous time step t — 1, but also the configuration at time 
stept-2 [109,3]. 

Nishinari et al. recently extended the BCA model, thereby allowing for slow-to- 
start effects with t> max > 1 cell/time step [110]. Their model contains a rule similar 
to the classic notion of slow-to-start rules, but now generalised for moving vehicles, 
leading to the terminology of a slow-to-accelerate rule. Taking the idea of antici- 
pation one step further, they also incorporated a driver's perspective, meaning that 
a vehicle will base its acceleration and braking decisions not only on the basis of 
its space gap and the anticipated speed of the vehicle ahead, but also on the space 
gap with the next leading vehicle (or even a vehicle located more downstream). As 
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a result, the model exhibits multiple metastable branches in the (k,q) diagram, as 
can be seen in Fig. 29. For the lowest metastable branch, vehicles inside jams will 
come to a complete stop. In contrast to this, vehicles will still be able to move for- 
ward inside jams for the higher branches. Note that depending on the strength of 
a local perturbation, traffic will shift from the highest branch to one of the lower 
branches. Finally, Nishinari et al. also combined the model with the classic STCA 
(see section 3.2.1), thereby allowing for stochasticity in both the acceleration and 
braking rules. 




Global density [vehicles^cells] 

Fig. 29. A (k,q) diagram of Nishinari et al.'s extended BCA model, with v max = 5 
cells/time step, AT = 1.3 s, AX = 7.5 m, and a driver's perspective of two vehicles 
ahead. The resulting diagram exhibits multiple metastable branches. Vehicles inside jams 
come to a complete stop only for the lowest metastable branch; for the higher branches, 
vehicles inside jams are still able to move forward. Depending on the strength of a local 
perturbation, traffic will shift from the highest branch to one of the lower branches (image 
reproduced after [110]). 



4 Multi-cell models 

Whereas all the previously discussed TCA models were based on a single-cell 
setup, this section introduces some of the existing multi-cell TCA models (still 
for single-lane traffic). In a multi-cell model, a vehicle is allowed to span a number 
of consecutive cells in the longitudinal direction, i.e., k > I cell. 

In the subsequent sections, we discuss several multi-cell TCA models encountered 
in literature. We first start with an overview of the artifacts that can be introduced 
when switching to a multi-cell setup. Subsequently, we describe three multi-cell 
TCA models, which have more intricate rule sets than the simple models of section 
3: 

• Helbing-Schreckenberg TCA (HS-TCA) 
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• Brake-light TCA (BL-TCA) 

• The model of Kerner, Klenov, and Wolf (KKW-TCA) 

Note that with respect to the measurements performed on the TCA models' lat- 
tices, we assume homogeneous traffic flows, i.e., all vehicles have the same length. 
This allows us, after suitable adjustment with the average vehicle length 7 = l u to 
express the global density as k g e [0, 1]. 



4. 1 Artifacts of a multi-cell setup 



It might seem that a translation of the classic STCA model (see section 3.2.1) into 
a multi-cell version would be straightforward. However, using a finer discretisation 
introduces a very specific artifact, i.e, hysteresis. In order to investigate this phe- 
nomenon, we have performed several experiments based on a multi-cell translation 
of the STCA model (now called the MC-STCA). In what follows, we assume a 
closed-loop lattice consisting of 10 5 cells. The simulations ran each for 5 x 10 5 
time steps, with AT = 1 s. 

Setting the slowdown probability to p = 0.5, the left part of Fig. 30 shows the result- 
ing (k,q) diagrams for different spatial discretisations, each time for homogeneous 
initial conditions. The average vehicle length was set to 7 G {2, 4, 8, 16, 32, 64} 
cells. In these experiments, we also scaled the maximum speed v max correspond- 
ingly (e.g., if 7 = 4 cells, then v max would become 5 x 4 = 20 cells/time step), 
as can be seen from the coinciding free-flow branches in the left part in Fig. 30. 
We also notice that an increase of the average vehicle length apparently results in 
a higher critical density, with an associated higher capacity flow. Furthermore, the 
flow seems to encounter a capacity drop at this critical dense. 
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Fig. 30. Left: several (k,q) diagrams of the MC-STCA, for I e {2, 4, 8, 16, 32, 64} cells 
and p = 0.5. As can be seen, an increase of the average vehicle length apparently results 
in a higher critical density, with an associated higher capacity flow (followed by a capacity 
drop). Right: the same setup for the MC-STCA, but now with a fixed 7 = 8 cells and 
Vmax = 5 x 8 = 40 cells/time step. The (k,q) diagrams depict the results of changing the 
slowdown probability p € {0.1,0.3,0.5,0.7,0.9}: an increase of p, leads to decrease of 
both the critical density and the capacity flow. 
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What causes this capacity drop ? To answer this question, we must first consider 
what happens in the deterministic case where p = 0. Here, our experiments have 
shown that there is no difference between a single-cell and a multi-cell setup. Set- 
ting p > 0, the randomisation rule R2, equation (47), introduces fluctuations in 
the high speeds of vehicles in free-flow traffic. However, these speed fluctuations 
are actually small compared to the vehicles' speeds themselves. Because of this 
limited influence, the free-flow branch of the (k,q) diagrams remains very stable. 
The smaller the discretisation, i.e., the larger the average vehicle length, the more 
stable the free-flow branch becomes for larger densities (note however that the ca- 
pacity drop gets less pronounced for increasing average vehicle lengths). This ca- 
pacity drop behaviour due to a stabilisation effect, is akin to the observations in the 
STCA's cruise-control limit (see section 3.2.2), and thus different from the VDR- 
TCA (see section 3.3.3), where a reduced outflow from a jam causes the drop in 
flow [19]. In contrast to this, random initial conditions or a superjam to start the 
simulations with, will always lead to the congested branch, thereby indicating a 
hysteretic phase transition. As the left part of Fig. 30 indicates, changing the dis- 
cretisation level of the STCA, by adjusting the average vehicle length and relatively 
keeping the same maximum speed, has only an effect on the length of the free-flow 
branch; the traffic dynamics in the congested regime remain the same. 

Holding 7 fixed at 8 cells and t> max = 5 x 8 = 40 cells/time step, the right part 
of Fig. 30 shows the resulting (k,q) diagrams for different values of the slowdown 
probability p e {0.1, 0.3, 0.5, 0.7, 0.9}. It is clear that an increase of p, leads to a 
decrease of both the critical density and the capacity flow. Note that the size of the 
capacity drop remains approximately the same for the different p. 

To conclude, we mention the work of Grabolus who performed extensive numeri- 
cal studies on the STCA. He also noted that it is possible to translate any multi-cell 
STCA variant into an equivalent single-cell STCA model, by suitably adjusting the 
values of the density and the maximum speed [100]. 

I Interestingly, the use of a smaller discretisation was already considered by 
Barrett et al. in the early course of the TRANSIMS project [111,1]. In their 
work, they introduce the terminology of multi-resolution TCA models, cor- 
responding to our multi-cell setup. Although they discuss several methods 
for integral refinements of the TCAs lattice, they do not make any mention 
of the observed hysteresis phenomenon introduced by a finer discretisation. 



4.2 Advanced multi-cell models 

Having discussed the repercussions of switching to a multi-cell setup, we now il- 
lustrate three TCA models that have more complex rule sets. We discuss their prop- 
erties by means of time- space diagrams, fundamental diagrams of global and local 
measurements, and histograms of the distributions of the space and time gaps. 
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4. 2. 1 The model ofHelbing and Schreckenberg ( HS-TCA ) 

Similar in spirit as the STCA (see section 3.2.1) and the ER-TCA (see section 
3.2.5), Helbing and Schreckenberg proposed their HS-TCA model in analogy with 
the optimal velocity model [112]. In fact, their model can be seen as a direct dis- 
cretisation of the OVM, with the following rule set: 

Rl: acceleration and braking 

Vi (t) <- Vi (t -l) + [a (V(g Si (t - 1)) - Vi (t - 1))J , (64) 
R2: randomisation 

£(t) <p^ Vi (t) <- max{0,Ui(*) - 1}, (65) 
R3: vehicle movement 

Xi(t)<-Xi(t-1)+Vi(t). (66) 

The function V(g Si ) in rule Rl, equation (64), is the discrete version of the optimal 
velocity function; it is specified in the form of a lookup table, containing speed en- 
tries for each space gap (see Table 1) and has the following meaning: higher values 
for the parameter a indicate an almost instantaneous adaptation of the vehicle's 
speed to the OVF, whereas lower values denote an increasing inertia and longer 
adaptation times [112]. However, as stated by Chowdhury et al. and Knospe et al., 
the role of a is a bit unclear as it does not exactly correspond to the timescale of the 
adaptation to the OVF (which is the case for the original optimal velocity model) 
[3,19]. Furthermore, certain values for a can, in combination with the OVF, lead 
to collisions between vehicles (because a reduces a vehicle's braking capability). 
Knospe et al. later provided the necessary conditions that guarantee collision-free 
driving, and avoid the possible backward moving of vehicles [19]. Note that, similar 
to the Fukui-Ishibashi models (see sections 3.1.2 and 3.2.3), vehicles are allowed to 
accelerate instantaneously in the HS-TCA model. The model is stochastic, in that 
it introduces randomisation by means of rule R2, equation (65), with £(£) e [0, 1[ 
a random number drawn from a uniform distribution. 

In Fig. 31, we have given two time-space diagrams of the HS-TCA for global den- 
sities k = 0.25 and k = 0.40 vehicles/cell. The length of a vehicle was / = 2 cells, 
p = 0.001, a = 1-r 1.3, v max = 15 cells/time step, AT = 1 s, and AX = 2.5 m. 
Due the small slowdown probability, the system dynamics are strongly determinis- 
tic, totally dependent on the initial (homogeneous) conditions. In the left diagram 
we can observe how vehicles can accelerate instantaneously when exiting a jam. 
Note that for higher densities, all jams become dense and compact, always con- 
taining stopped vehicles, as is depicted in the right diagram. Because of the non- 
linearity introduced by the discretised optimal velocity function, all tempo-spatial 
patterns in the system are of a chaotic nature (i.e., nonlinear with stochastic noise) 
[19]. 

The (k,v s ) and (k,q) diagrams in Fig. 32 are based on local and global measure- 
ments. A feature of these diagrams is that the local measurements tend to form 
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Table 1 

A possible optimal velocity function (OVF) for the TCA model of Helbing and Schreck- 
enberg (HS-TCA). The OVF is represented as a table, giving the optimal speed V(g Si ) 
(expressed as cells/time step) associated with each possible space gap g Si (expressed as a 
number of cells). 




Fig. 31. Typical time-space diagrams of the HS-TCA model, with I = 2 cells, p = 0.001 
a = 1-7- 1.3, and v max = 15 cells/time step. The shown closed-loop lattices each contain 
300 x2 = 600 cells, with a visible period of 580 time steps. The global density k was set to 
0.25 vehicles/cell (left) and 0.40 vehicles/cell (right). The formation of congestion waves 
leads to dense, compact jams containing stopped vehicles. Vehicles strive to decelerate 
smoothly, but are allowed to accelerate instantaneously when exiting jams fronts. 

clusters around certain space-mean speeds (see the left part of Fig. 32): these clus- 
ters correspond to the speeds dictated by the discretised optimal velocity function 
of Table 1, each time associated with an average space gap corresponding to the in- 
verse of the locally measured density. As a result, the (k,q) diagram in the right part 
of Fig. 32 shows several branches, each one with a different OVF speed. The lowest 
branch corresponds to the speed of the backward propagating waves, i.e., the jam 
speed. Even more striking, is that from a certain finite density k <^ 1 vehicle/cell 
on, all vehicles always come to a full stop and the flow in the system becomes zero 
[112]. 

To conclude our discussion of the HS-TCA, we give the histograms of the distribu- 
tions of the space and time gaps in the left and right parts, respectively, of Fig. 33. 
The most prominent features of these histograms, are that (i) there exist small clus- 
ters of probability mass between certain space gaps (i.e., 15 - 20, 25 - 25, and 35 
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Fig. 32. The (k,v s ) (left) and (k,q) (right) diagrams for the HS-TCA, obtained by local 
and global measurements. The local measurements tend to form clusters around certain 
space-mean speeds, corresponding to the speeds dictated by the discretised optimal velocity 
function of Table 1 . These clusters are visible in the right diagram as branches with different 
slopes. Remarkably, from a certain finite density k <C 1 vehicle/cell on, all vehicles always 
come to a full stop and the flow in the system becomes zero. 



- 40 cells), corresponding to groups of vehicles, (ii) for higher densities, we can 
observe a spread-out cluster of probability mass, corresponding to the lowest local 
measurements in the left part of Fig. 32, and (iii) in contrast to the previous TCA 
models, the median of the time gap for the HS-TCA is already very small for den- 
sities k < 0. 1 . 





Fig. 33. Histograms of the distributions of the vehicles' space gaps g s (left) and time gaps 
gt (right), as a function of the global density k in the HS-TCA. The thick solid lines denote 
the mean space gap and median time gap, whereas the thin solid line shows the former's 
standard deviation. The grey regions denote the probability densities. 
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The HS-TCA might seem an interesting improvement, as its being based on 
a discretisation of the optimal velocity model. But although its authors state 
that it "reproduces many of empirically observed features" [112], Knospe 
et al. showed several shortcomings in the model [19]: care must be taken 
to avoid collisions, and the model fails to reproduce the synchronised- flow 
regime entirely. This latter can be understood by looking at the dense, com- 
pact structure of jams in the time-space diagrams of Fig. 31, and the occur- 
rence of branches with distinct speeds as in the right part of Fig. 32. 



4.2.2 Brake-light TCA (BL-TCA) 

Recently, an interesting idea was pursued by Knospe et al.; their TCA model in- 
cludes anticipation effects, introduced by equipping the vehicles with brake lights 
[113]. The focus of this (and the following) TCA model lies in a correct repro- 
duction of the three phases of traffic as introduced by Kerner et al. [58,1]. In a 
sense, the BL-TCA incorporates many of the features encountered in previously 
discussed single-cell TCA models. First of all, the BL-TCA has randomisation for 
spontaneous braking. Secondly, it has slow-to- start behaviour for the capacity drop 
and hysteresis phenomena. Moreover, it incorporates anticipation which can lead 
to a stabilisation of the free-flow branch. Finally, it includes elements for reproduc- 
ing synchronised traffic. These latter two aspects clearly go beyond the standard 
incentive if drivers to avoid collisions. As such, it is the desire for smooth and 
comfortable driving (which resembles human behaviour), is responsible for the oc- 
currence of traffic states like e.g., synchronised traffic [40]. To achieve all this, the 
rule set of the BL-TCA becomes quite complex, in comparison with some of the 
more standard single-cell TCA models of section 3: 

R0: determine stochastic noise 



6 i+1 (t-l) = l A t hi (t- 1) <t Si (t- 1) 
=>• p(t) <- p b , 



Vi(t - 1) = 

=>• p(t) <- p , 

else 

=>• P(t) <- Pd, 
bi(t)<-0, 
Rl: acceleration 



(67) 



(6i(t - 1) = A 6 <+1 (t-l)=0) 
v t hi (t)>t Si (t) 

Vi(t) <- mm{Vi(t + 1), Umax}, 



(68) 
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R2a: determine effective space gap 



( 69 ) 

max{min{^ +1 (t - l),g 8i+l (t - 1)} -# Ssecurity , 0}, 



anticipated speed of leading vehicle 



R2b: braking 



Vi(t)^mm{vi(t),gl(t)}, 

Vi(t)<Vi(t-l) (70) 
=*► 6i(*) <- 1, 



R3: randomisation 



(71) 



£(f) < p(f) =► 
p(t)=p b A = Vi(t - 1) + 1 

=► hit) <- l, 

<- max{0,fi(t) - 1}, 

R4: vehicle movement 

Xi (t) ^ Xiit-^+v^t), (72) 

where fej(i) denotes the state (0 or 1) of the brake light of the i th vehicle at time step 
t, thi = 9si/ v i an d t Si = min{vi, h} with h the interaction range of the brake light. 
As such, ^ s m ^ time to reach the leading vehicle, which gets compared with an 
interaction horizon t Si that depends on the speed vi and is constrained by h. If the 
leading vehicle is far away, its brake light should not influence the following vehi- 
cle. Furthermore, rule R0 also takes into account that drivers are more alert when 
they are travelling at high speeds. The slowdown probability p in rule R0, equation 
(67), corresponds to either the braking probability p b , the slow-to-start probability 
Po, or the classic slowdown probability p d for decelerations. Finally, g* (t) in rules 
R2a and R2b, equations (69) and (70), respectively, denotes the effective space gap, 
based on the anticipated speed of the leading vehicle and taking into account a se- 
curity constraint g Sssxmit • Just as the previous TCA models, the BL-TCA is stochas- 
tic, in that it introduces randomisation by means of rule R3, equation (71), with 
£(£) G [0, 1[ a random number drawn from a uniform distribution. If a vehicle was 
in the process of braking due to the previous rules, then its brake light 6, is turned 
on. Note that Knospe also extended the BL-TCA with rules that allow asymmetric 
lane changing on a two-lane road (unidirectional), incorporating a right-lane pref- 
erence as well as an overtaking prohibition on the right lane. As such, the model 
correctly reflects the density inversion phenomenon (see also section 5.1) [40,1 14]. 
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In the remainder of this discussion, we set p b = 0.94, p = 0.5, pd = 0.1, h =6 
time steps, g Sxciu ^ = 7 cells, t> max = 20 cells/time step, with a vehicle length of 
/ = 5 cells, AT = 1 s, and AX = 1.5 m [113,19]. With respect to the calibration 
of the BL-TCA model's parameters, Knospe et al. provide a nice overview, giving 
intuitive analogies for each of these parameters (e.g., p is associated with the speed 
of the backward propagating waves) [19]. 

In Fig. 34, we have given two time-space diagrams of the BL-TCA for global densi- 
ties k = 0.25 and k = 0.40 vehicles/cell. As can be seen in the time-space diagram 
in the left part, the anticipation and synchronisation phenomena lead to forward 
propagating density waves, where vehicles carry the density downstream. Going to 
higher densities, we can see stable jams, indicative of the wide-moving jam phase 
(see also Kerner's three-phase traffic theory [58,1]). 




Fig. 34. Typical time-space diagrams of the BL-TCA model (refer to the text for the used 
parameter values). The shown closed-loop lattices each contain 300 x5 = 1500 cells, with 
a visible period of 580 time steps. The global density k was set to 0.25 vehicles/cell (left) 
and 0.40 vehicles/cell (right). The visible forward propagating density waves are a result 
of the anticipation and synchronisation phenomena. At higher densities, stable jams occur, 
indicative of the wide-moving jam phase. 

Looking at the (k,v s ) and (k,q) diagrams in Fig. 35, we can use the local measure- 
ments to discriminate between the free-flow (o), synchronised-flow (■), and jammed 
regimes (*). The synchronised regime is visible as a wide scatter in the data points, 
having various speeds but relatively high flows. The data points in the wide-moving 
jam correspond to Kerner's so-called line J [58,1]. The use of a finer discretisation 
can lead to metastable states (see section 4.1), but as Knospe et al. note, the slow- 
to-start behaviour in rule R0, equation (67), is necessary in order to produce the 
correct speed of the backward propagating wave, as a result of a reduced outflow 
from a jam [19]. 

Finally, Fig. 36 depicts the histograms of the distributions of the space and time 
gaps in the left and right parts, respectively. In contrast to the HS-TCA, there are 
no more clusters for the space gap (see left part of Fig. 33), but rather a smooth re- 
gion of probability mass: as the global density of the system increases, the average 
space gap diminishes continuously and monotonically. The observations for the dis- 
tributions of the time gaps correspond to those encountered in literature [113,19]: 
from the right part of Fig. 36, we can see a wide range of probability mass at low 
densities (free-flow traffic), corresponding to a wide distribution of time gaps. At 
intermediate densities (synchronised flow), the distribution tends to peak, leading 
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Fig. 35. The (k,v s ) (left) and (k,q) (right) diagrams for the BL-TCA model, obtained by 
local and global measurements. The local measurements discriminate between the free-flow 
(o), synchronised-flow (•), and jammed regimes (*). The synchronised regime is visible as 
a wide scatter in the data points, having various speeds but relatively high flows. The data 
points in the wide-moving jam correspond to Kerner's line J. 

to a small dense cluster at approximately k =0.15 vehicles/cell, with a median time 
gap of 1 time step. Finally, at higher densities (jammed traffic), the distribution of 
the time gaps gets more peaked, as is illustrated by the narrowing of the grey region 
of probability mass. 




Fig. 36. Histograms of the distributions of the vehicles' space gaps g s (left) and time gaps 
gt (right), as a function of the global density k in the BL-TCA model. The thick solid lines 
denote the mean space gap and median time gap, whereas the thin solid line shows the 
former's standard deviation. The grey regions denote the probability densities. 

4.2.3 The model ofKerner, Klenov, and Wolf(KKW-TCA ) 

Based upon the BL-TCA of Knospe et al., Kerner, Klenov, and Wolf (KKW) re- 
fined this approach by extending it. Their work resulted in a family of models that 
incorporate the notion of a synchronisation distance for individual vehicles [115]. 
Derived from this model class, Kerner et al. proposed discretised versions in the 
form of traffic cellular automata models. In this report, we consider the KKW-1 
TCA model, of which the complex rule set is as follows [116]: 
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Rla: determine synchronisation distance 

D .(t)^D Q + D l v i (t-l), 
Rib: determine acceleration and deceleration 

Vi(t - 1) < v i+ i(t - 1) =>- A acCi (t) <- a, 
Ui(t - 1) = - 1) =^ A acCi (t) <- 0, 

Ui(t - 1) > u i+ i(t - 1) =>• A acCi (t) < 6, 

Rlc: determine desired speed 

g 9i (t-l)>(D i (t)-l i ) 

=>• v deSi (t) <- - 1) + a, 

fc(«-l)<(A(t)-li) 

=>• ^des, (*) <- Ui (* - 1) + A acCi (t) , 

Rid: determine deterministic speed 

Vi(t) <- max{0,min{t; max ,5( Si (t),t; d e Sl (t)}}, 
R2a: determine acceleration probability 



Vi(t) < v p =>p a (t) ^p ai , 

Vi(t) > V p =>p a (t) ^Pa 2 , 

R2b: determine braking probability 

Vi(t) = 0^p b (t) <-p , 
Vi(t) >0^p b (t) <-p d , 
R2c: determine stochastic noise 



£(*) < Pa(*) 



Pa(0 < < Pa(t)+Pb(t) =}► 7fc(*) 



>Pa(*)+P&(*) 

R2d: determine stochastic speed 



rh(t)<-0, 



Vi(t) <- max{0, min{w max , ^(t) + r]i(t),Vi{t) + a}}, 
R3: vehicle movement 



Xi(t) <- - 1) + ^(t). 
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As can be seen from this overview, the KKW-TCA model's rule set is mainly com- 
posed of a deterministic part (rules Rla - Rid) and a stochastic part (rules R2a - 
R2d). In the deterministic part, the synchronisation distance D,i is computed first 
with rule Rla, which uses a linear function (other forms, e.g., quadratic functions, 
are also possible). The parameters Dq and D\ need to be estimated. Rule Rlc de- 
termines the desired speed v^:. the first part of the rule allows the vehicle to ac- 
celerate, whereas the second part of the rule uses an acceleration A acCi defined by 
rule Rib (a and b are parameters denoting the acceleration, and respectively brak- 
ing, capabilities). As such, a vehicle will tend to adapt its speed to that of its direct 
frontal leader, whenever the vehicle is within a zone of interaction (i.e., the syn- 
chronisation distance). The deterministic speed is then computed by means of rule 
Rid, which takes into account the maximum speed v max , the space gap g Si to avoid 
a collision, and the previously computed desired speed of rule Rlc. 

In the stochastic part for computing the speed, a randomisation is introduced in rule 
R2d by means of a stochastic acceleration r\i. The values of r\i are obtained in rule 
R2c with probability p a for accelerating, and probability p b for braking. The former 
is dependent on the vehicles computed deterministic speed and the parameters v p , 
p ai , and p a2 with p ai > p ai and p a , + Pa 2 < 1 • The latter, p b is dependent on the 
vehicles computed deterministic speed and the slowdown probability pd and the 
slow-to-start probability p withp > Pd- 

In the remainder of this discussion, we set D = 60, D\ = 2.55, a = b = 1, 
v p = 28, p ai = 0.2, p a2 = 0.052, p = 0.425, pd = 0.04, v max = 60 cells/time step, 
with a vehicle length of / = 15 cells, AT = 1 s, and AX = 0.5 m [19]. 

Considering the KKW-TCA models' time-space diagrams in Fig. 37, we can see 
that, in contrast to the BL-TCA (see section 4.2.2), there are less spontaneous for- 
mations of small traffic jams. The forward propagating density waves in Fig. 34 are 
absent in the KKW-TCA model. However, the two models show good correspon- 
dence with respect to the speed of the backward propagating waves. 




Fig. 37. Typical time-space diagrams of the KKW-TCA model (refer to the text for the used 
parameter values). The shown closed-loop lattices each contain 300 x 15 = 4500 cells, with 
a visible period of 580 time steps. The global density k was set to 0.25 vehicles/cell (left) 
and 0.40 vehicles/cell (right). Note the stable flow of vehicles surrounding the dense and 
compact superjams. 

Similar as in the BL-TCA model's effective space gap g*.(t), the synchronisation 
distance D is responsible for producing the typical two-dimensional scatter in the 
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(k,v s ) and (k,q) diagrams in Fig. 38. When a driver who is within the synchroni- 
sation distance adapts the vehicle's speed, the only factors taken into account are 
the current speed of the direct frontal leader and a safety criterion (in the form of 
the current space gap); it is this effect that produces the scatter in the data, because 
the exact specification of this speed is absent. In both diagrams of Fig. 38, the local 
measurements discriminate between the free-flow (o), synchronised-flow (■), and 
jammed regimes (*). One of the major differences between these two models, is 
that the flow in the synchronised regime is almost a factor two larger for the KKW- 
TCA than the BL-TCA. The KKW-TCA also experiences a capacity drop similar 
as in the BL-TCA, but also undergoes an abrupt transition when going from the 
synchronised-flow to the wide-moving jam regime around a global density of some 
0.4 vehicles/cell (see the left part of Fig. 38). Because the model is built around the 
assumption that vehicles tend to approximate the behaviour of their direct leader 
within a certain synchronisation distance, the resulting traffic regimes correspond 
well to Kerner's empirical observations [58,1]. 




Fig. 38. The (k,v s ) {left) and (k,q) {right) diagrams for the KKW-TCA model, obtained by 
local and global measurements. The local measurements discriminate between the free-flow 
(o), synchronised-flow (•), and jammed regimes (*). The synchronised regime is visible as a 
wide scatter in the data points, having various speeds but flows comparable to the capacity 
flow. 

In Fig. 36, we have depicted the histograms of the distributions of the space and 
time gaps in the left and right parts, respectively. The distributions are similar to 
those of the BL-TCA, but there are some important differences. With respect to 
the space gaps in the left part of Fig. 36, there is a high variance in the jammed 
regime, due to the fact that there are vehicles in free-flow traffic, as well as inside 
the wide-moving jams (although most of the probability mass is assigned to the 
zero space gap inside the dense jams). Considering the time gaps in the right part 
of Fig. 36, we can see that they always form a tight cluster around the median of the 
distribution, indicating very narrow distributions with an pronounced peak. This is 
completely different behaviour than in the BL-TCA model (see the right part of 
Fig. 36). The main reason is probably due to the lack of an anticipation effect in 
the KKW-TCA model. Even more severe, is the fact that the KKW-TCA model, 
despite its elaborate construction based on a synchronisation distance, completely 
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fails to describe the microscopic structure of motorway traffic. The BL-TCA model 
however succeeds in having a good fit on both macroscopic and macroscopic scales, 
as stated according to Knospe et al. [114,19]. 




Fig. 39. Histograms of the distributions of the vehicles' space gaps g s (left) and time gaps 
g t (right), as a function of the global density k in the KKW-TCA model. The thick solid 
lines denote the mean space gap and median time gap, whereas the thin solid line shows 
the former's standard deviation. The grey regions denote the probability densities. 

5 Multi-lane traffic, city traffic, and analytical results 

In this final section on traffic cellular automata models, we take a look at some other 
aspects related to TCA models. We first discuss some properties and methodologies 
for modelling multi-lane traffic in the context of a cellular automaton, after which 
we briefly consider several approaches for dealing with city traffic. The final part of 
the section concludes with an overview of different analytical treatments of TCA 
models. 

5. 1 Multi-lane traffic 

In this section, we briefly discuss some properties and methodologies for modelling 
multi-lane traffic in the context of a cellular automaton. To this end, we illustrate 
the types of lane changes that are possible, then discuss the general setup for a 
lane-changing model. We conclude with a short overview on the implementation 
of lane-change rules and explain the phenomenon of ping-pong traffic, an artifact 
introduced by an inferior implementation. 

5.1.1 Types of lane changes 

In general, there are two types of lane changes identified: mandatory lane changes 
(MLC) and discretionary lane changes (DLC) [1]. In the former case, a vehicle is 
obliged to execute a lane change, e.g., because it needs to exit the motorway at an 
off-ramp, or because the vehicle is by law obliged to drive in the right shoulder 
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lane. In the latter case, a vehicle changes a lane at its own discretion, e.g., when 
approaching and overtaking a slow-moving leading vehicle. 

With respect to the rules for lane changing, there are also two approaches: symmet- 
ric and asymmetric. In the US, the symmetric approach is more applicable: this is 
embodied by the fact that motorways have a large number of lanes (i.e., more than 
three), with vehicles driving at lower speeds (e.g., 60 miles/hour, corresponding to 
some 100 kilometres/hour), effectively using all lanes more homogeneously. Such 
a system is typically called "keep-your-lane", as frequent lane changes are discour- 
aged. In contrast to this, people in most European countries are obliged by law to 
drive on the outer right shoulder lane whenever possible. Motorways have fewer 
lanes (typically either two or three, unidirectional), operating at higher speeds of 
e.g., 120 kilometres/hour. In addition, most of these countries have instituted an 
overtaking prohibition on the right lane, with large trucks restricted to the two most 
right lanes. 

With respect to this latter system of asymmetric lane changes, the phenomenon 
of density or lane inversion plays an important role, especially on the numerous 
2x2 motorways in Europe (see also the beginning of section 3 for a discussion of 
this phenomenon). Another aspect that has a significant influence, is the change of 
driver behaviour, e.g., near on-ramps. Here, drivers might avoid the shoulder lane 
to allow traffic to enter, or because of their increased attention, they might induce 
a more subtle effect such as the capacity funnel (see also our discussion in [51] for 
more details on this phenomenon). 

5.1.2 General setup for lane changing 

Deciding on whether or not to perform a lane change, is typically split in two sep- 
arate steps: first, a vehicle checks if it is desirable to change lanes, i.e., making the 
distinction between a mandatory or discretionary lane change. If a lane change is 
indeed desirable, then the second step proceeds to check whether or not such a lane 
change can be performed at all with respect to safety and collision avoidance. Thus, 
there is a check for gap acceptance. 

One of the first approaches to model such lane-changing behaviour an a two-lane 
road in a TCA model, is due to Nagatani. His work was based on the deterministic 
CA-184 model (see section 3.1.1) [117]. One of the artifacts of his lane-changing 
rules, was the existence of states in which blocks of vehicles alternated from one 
lane to another, without moving at all. To circumvent this problem, Nagatani ran- 
domised the lane-changing behaviour [118]. Rickert et al. later applied this lane- 
changing methodology, by extending the STCA model (see section 3.2.1) to handle 
two-lane unidirectional traffic [1 19]. Wagner et al. later assessed the previous work 
of Rickert et al., concluding that it did not capture certain aspects (e.g., density in- 
version) of traffic flows very well [120]. To this end, they built upon the previous 
work, adding a more specialised security constraint that takes into account the fact 
that vehicles should also consider the following vehicles in the target lane, thereby 
avoiding severe disruptions. As a final comment, they state that the lane-changing 
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rules in a TCA model typically do not provide a realistic microscopic model, but 
they rather lead to a good correspondence with respect to observed macroscopic 
features (e.g., the frequency of lane changes). 

In order to address the correct reproduction of the density inversion phenomenon, 
Nagel et al. artificially introduced a slack parameter, capturing the inclination of a 
driver to change back to the right lane. They furthermore also provided an extensive 
classification of some 10 lane changing rules and criteria encountered in literature 
[56]. Another excellent overview of multi-lane traffic is given by Chowdhury et al. 
[3]. 

As all the previous work dealt with unidirectional roads, it seems logical to con- 
sider bidirectional traffic, i.e., traffic with adjacent but opposing lanes. Simon and 
Gutowitz were among the first to consider a TCA model of such traffic, with vehi- 
cles driving on two lanes [121]. Central to their approach, is the notion of a local 
density that each driver must assess before attempting to complete an overtaking 
manoeuvre. When a driver encounters a slower moving vehicle, a check is made 
whether or not there is enough space in front of this leading vehicle (this is the lo- 
cal density). If the check is positive, then a lane change can be performed (under the 
condition of course that there is a safe gap in the opposing lane). With this scheme 
in mind, high density traffic thus excludes such overtaking manoeuvres, due to the 
fact that the local density is too low to complete them. 




Note that some authors, e.g., Gundaliya et al. [122], Mallikarjuna and Ra- 
machandra Rao [123], use a peculiar variant of a multi-lane setup. Their 
models have essentially a multi-cell structure, but now the multi-cell con- 
cept is extended in the lateral direction. So cells not only get smaller, but 
also 'thinner', allowing variable -width vehicles, e.g., motor cycles that can 
more easily pass other vehicles in the same lane. In our opinion, this leads to 
unnecessary complexity, giving little benefits. In fact, we believe that such a 
scheme directly opposes the idea behind a CA model, as explained at the in- 
troduction of this report. We strongly feel that heterogeneity in a TCA model 
should only be incorporated by means of different lengths, maximum speeds, 
acceleration characteristics, anticipation levels, and stochastic noise for dis- 
tinct classes of vehicles and/or drivers. Any other approach would be better 
off with a continuous microscopic model. 



5.1.3 Implementation of lane-changing rules and the phenomenon of ping-pong 
traffic 

The basic implementation of a lane-changing model in a TCA setting, leads to two 
substeps that are consecutively executed at each time step of the CA: 

• first, the lane-changing model is executed, exchanging vehicles between laterally 
adjacent lanes, 

• then, all vehicles are moved forward (i.e., longitudinal) by applying the car- 
following part of the TCA model's rules. 
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One immediate result from this approach, is that a lane change in a TCA model 
is completed within one time step (i.e., AT). This is in contrast to real-life traffic, 
where lane changes have a duration of several seconds [56]. 

For more than two lanes, care must be taken to avoid so-called scheduling conflicts 
during the first substep. Consider for example three lanes, with two vehicles driving 
in the outer left, respectively outer right, lane at the same longitudinal position. If 
the cell in the middle lane is empty, then the vehicles may decide to move to this lo- 
cation, resulting in a lateral collision. In order to compensate this, one possibility is 
to choose a vehicle at random (or by preference), thereby allowing it to perform its 
requested lane change. Another possibility is to perform left-to-right lane changes 
in even time steps, and right-to-left lane changes in odd time steps. 

As hinted earlier, the 'correctness' of a lane-change model should be judged on the 
basis of certain macroscopic observations. Examples of these are the frequency of 
lane changes with respect to different densities, the capacity flows for all lanes sep- 
arately and combined, the critical density at which a breakdown occurs in each of 
the lanes, . . . Good indicators can be found in the many small fluctuations typically 
exhibited by multi-lane TCA models, instead of the large jams in single-lane traf- 
fic. Traffic flows get more fluid if vehicles are allowed to pass moving bottlenecks 
[120,56]. However, under certain conditions, Helbing and Huberman have shown 
the existence of coherent states, where vehicles' speeds are synchronised across ad- 
jacent lanes. For heterogeneous traffic flows, this can lead to a moving 'solid block' 
of vehicles [124]. 

When implementing lane-change rules in a TCA model, care must however be 
taken that the implementation does not introduce any unrealistic artifacts. A promi- 
nent example of this, plaguing many TCA models, is a phenomenon called ping- 
pong traffic. Nagatani was among the first to observe this peculiar behaviour of ve- 
hicles in traffic flows (see section 5.1.2). In ping-pong traffic, vehicles typically al- 
ternate between lanes during successive time steps. As explained earlier, one way to 
resolve this behaviour is by randomising the lane-change decision, thereby quickly 
destroying any such artificial patterns [118,119]. 

5.2 City traffic and intersection modelling 

When modelling city traffic, essentially two approaches can be followed: either 
the entire road network is considered as a two-dimensional lattice (i.e., a grid), or 
each road in the network is a single longitudinal lattice (single- or multi-lane) with 
explicitly modelled intersections. The former was historically used in the context of 
phase transitions in a CA, whereas the latter is more applicable to describe real-life 
traffic flows in populated cities. 

In this section, we illustrate both approaches, starting with a classic grid layout as 
embodied by the Biham-Middleton-Levin (BML) and Chowdhury-Schadschneider 
(ChSch) TCA models, after which we briefly comment on explicit descriptions of 
intersections in TCA models. 
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5.2.1 Grid traffic 



The first model of 'city traffic' was proposed by Biham, Middleton, and Levine 
(BML). It was developed around the same time Nagel and Schreckenberg presented 
their STCA (see section 3.2.1). The BML-TCA, is a two-dimensional model that 
describes traffic on a square grid in a toroidal setup (i.e., opposing sides are iden- 
tified), with vehicles distributed randomly over the lattice [125]. The model is in 
fact a very simplistic model, in that it assumes that all vehicles either move from 
the south to the north direction, or from the west to the east. Each cell of the lattice 
is assumed to contain a traffic light, in the sense that all west-east vehicles try to 
move during even time steps, and all south-north vehicles during odd time steps 
(thus f max = 1 cell/time step for all vehicles). The BML-TCA constitutes a fully 
deterministic model, where the only randomness is introduced through the initial 
conditions. Note that its one-dimensional version corresponds to the CA-184 and 
the TASEP (see sections 3.1.1 and 3.2.4). 

Depending on the global density of vehicles in the lattice, the model results in 
two distinct traffic regimes, with a sharp first-order phase transition between them. 
The first regime, i.e., free-flow traffic, corresponds to a state with alternate moving 
vehicles (i.e., west-east and south-north moving); an example is depicted in the left 
part of Fig. 40. In the congested regime, a self-organised global cluster emerges, 
completely composed of blocked vehicles (see e.g., the right part of Fig. 40). When 
the phase transition between both regimes occurs, the space-mean speed changes 
abruptly from one to zero cells/time step [125,126]. Fukui and Ishibashi studied 
the repercussions of a local disruption in the lattice (e.g., a crashed vehicle that 
remains stopped for an eternal period), and found that it provides the seed of a 
growing global cluster [127]. Biham et al. also considered a less restrictive version 
of the above model, in which now all vehicles try to move at each time step. In case 
of conflicts between a west-east and a south-north vehicle, one of them is chosen 
at random. Another variation considers also opposing traffic, which can lead to 
gridlocked situations where no vehicles are able to move at all. A generalisation of 
the BML-TCA, was provided by Freund and Poschel who consider a similar setup, 
but now with traffic moving in all four directions [128]. Finally, Shi is able to obtain 
analytical expressions for the critical densities at which the previously mentioned 
phase transitions occur [129]. 

In the work of Chowdhury et al., a comprehensive overview is given, describing 
extensions to the BML-framework [3]. This overview includes asymmetric distri- 
butions of the west-east and south-north vehicles, unequal maximum speeds, two- 
level crossings (where two vehicles can share the same cell), faulty traffic lights 
(here, either a west-east or south-north vehicle is chosen at random to occupy a 
cell, irrespective of the current time step), road blocks, line- and point-defects (i.e., 
a crowded 'street' of the model, corresponding to a dense horizontal or vertical row 
of cells), random turning of vehicles, cut-off streets (similar to a row of two-level 
crossings), and so forth and so on. 

Chowdhury and Schadschneider later extended the BML-TCA model to incorpo- 
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Fig. 40. Left: snapshot of the spatial structure in the BML-TCA for k = 0.25. In this 
free-flow regime, all vehicles move alternatingly, with the right-oriented arrows denoting 
west-east travelling vehicles, and the upward-oriented arrows denoting south-north trav- 
elling vehicles. Right: same setup as before, but now for k ss 0.4082. In this congested 
regime, a global cluster emerges, completely composed of blocked vehicles. Bottom: an 
overview of the ChSch-TCA, showing the street segments of finite length between the 
BML-TCA's original intersections. The first two images are reproduced after [125], the 
third after [98]. 

rate randomisation effects like in the STCA model, having the result that jamming 
can now occur spontaneously [130]. Their model furthermore contains street seg- 
ments of finite length between the cells, with vehicles driving according to the 
STCA's rules on these streets. The original cells in the BML-TCA model form the 
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signallised intersections of the Chowdhury-Schadschneider model (ChSch-TCA), 
as can be seen in the bottom part of Fig. 40. At sufficiently large densities, a transi- 
tion can occur that leads to a self-organising state of completely gridlocked traffic. 
Barlovic later provided a solution to this problem, making the model well-suited 
for assessing the results of different traffic light control policies in a city [98]. 



5.2.2 Explicit intersection modelling 

In contrast to the previous section were all traffic operations were essentially de- 
fined on a two-dimensional lattice, it is also possible to consider a complete road 
network, consisting of separate links that are connected to each other by means of 
intersections. These intersections can either be signallised, or unsignallised, turn- 
ing priorities can be defined, as well as different geometrical layouts (e.g., round- 
abouts). 



Road networks based on the above assumptions, typically combine a set of basic 
building blocks. As such, the network is logically decomposed in a set of nodes 
and links. The former denote the intersections, whereas the latter can, depending 
on the implementation, refer to individual lanes, a group of adjacent lanes, or even 
a road with two-way traffic. In general, traffic operations on motorways are primar- 
ily influenced by the behaviour of vehicles on links, i.e., their car-following and 
lane-changing behaviour. Conversely, traffic operations in cities and denser street 
networks, are primarily defined by the behaviour of vehicles at intersections, i.e., 
queueing delays at traffic lights, priority turns, ... In many cases, the intersection 
logic is simplified, such that all decisions (conflict resolving et cetera) are taken 
before a vehicle enters the intersection [14]. 

Several non-exhaustive examples include the work of Esser and Schreckenberg 
with applications to the city of Duisburg [131], the work of Simon and Nagel who 
primarily focussed on single-lane traffic in combination with several setups for con- 
trolling traffic lights, applying their work to the city of Dallas (different links have 
different slowdown probabilities associated with them, thus enabling to model dif- 
ferent street capacities) [132], the work of Diedrich et al. who consider the effects 
of various implementations of on- and off -ramps in the classic STCA model [133], 
and all the references on TRANSIMS, the travel behaviour in Switzerland, the re- 
gion of Dallas, the city of Portland, and the city of Geneva (where all intersections 
are replaced by generalised roundabouts), mentioned in our discussion in [1]. 
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All these examples have in common that they are based on simple building 
blocks. Despite this elegance, most of them however, do not provide satisfac- 
tory information regarding the calibration and validation of their underlying 
models (this for example with respect to the correct observed queueing de- 
lays at intersections). A popular technique is to use sources and sinks, where 
vehicles are added and removed, allowing tuning of the simulator in order 
to agree with incoming on-line measurements. Clearly, we feel that besides 
a need for elaborate descriptions of the employed models, there is perhaps 
even a bigger need for correct information with respect to these models' fi- 
delity and accuracy. 



5.3 Analytical results 

Because most studies based on TCA models heavily rely on numerical simulations, 
this creates the danger of introducing artifacts (e.g., finite-size effects) that obscure 
the true dynamics of the systems under consideration. Although most of these prob- 
lems should resolve in the so-called thermodynamic limit where Kc, T mp — > +00 
(i.e., a lattice with infinite length considered over an infinite time period), resort- 
ing to this approach is computationally not feasible. As a result, researchers have 
focussed on analytical methods. Except for the most trivial cases with a determin- 
istic (i.e., noiseless) TCA model, these analytical methods most of the time provide 
approximations at best. 

In this section, we illustrate several of these analytical methods encountered in lit- 
erature. Our discussion focusses on the concept of a mean-field theory, after which 
we elaborate on some of its improvements that lead to better agreement with nu- 
merical results. 

I Note that other avenues for analytical treatments of CA models, and TCA 
models in particular, are also explored. In this section, we will however not 
go into detail about them. For more information, we refer the reader to the 
interesting work of Fuks and Boccara [134,135,136,137,138]. 

5.3.1 Mean-field theory 

As mentioned in the introduction of this section, for the case of arbitrary v max and 
p = or p = 1, or for t> max = 1 cell/time step, the analytical solution of the re- 
sulting TCA model is exactly known. This solution, expressed as its (k,q) diagram, 
corresponds to the set of diagrams as depicted in Fig. 10 (see section 3.1.2) for the 
DFI-TCA. 

The problem is to find an analytical description of how the system evolves in time 
through the state space, i.e., what are the occurring configurations ? The evolution 
of a system, can be described by what is called a master equation. For cellular 
automata, this equation is a first-order differential equation, describing the change 
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in probability of a system's lattice to be in a certain configuration. The downside is 
that, in general, this master equation can not be solved exactly. 

For the TASEP model (see section 3.2.4) with open boundary conditions and ran- 
dom sequential update, the master equation can be solved exactly [139,18]. In a 
first step, the master equation is elegantly written in vector form, comprising a 
transfer matrix that contains the time-evolution of the probabilities. By assuming 
the matrix-product ansatz (MPA) formalism, the transfer matrix can be rewritten 
as a product of local transfer matrices, operating on sets of cells. This provides a 
algebra that can be solved exactly, thereby solving the TASEP analytically. Note 
that for the TASEP with a parallel update however, obtaining the exact solution is 
difficult, because no simple MPA decomposition into local matrices is possible. 

In contrast to this promising result, obtaining an analytical solution becomes harder 
to even intractable for the STCA model (see section 3.2.1) with v max > 1 cell/time 
step and < p < 1 . In the master equation, probabilities of cluster of cells will oc- 
cur, making its solution very hard [54]. One well-known method that is suitable for 
dealing with many-particle systems in statistical mechanics, is the construction of a 
mean-field theory (MFT) of the model. Such a MFT can provide an approximation 
of the master equation; in some cases, the MFT turns out to be an exact solution. 

The idea behind a MFT, is that all correlations between neighbouring cells are 
neglected. For TCA models, such a site-oriented mean-field theory (SOMF) as- 
sumes that all cluster probabilities are replaced by single cell probabilities. The 
MFT now replaces the effects of these individual cells with an average effect (the 
'mean field'), which simplifies computations considerably. When translating the 
STCA's rules Rl - R3, i.e., equations (46) - (48), Rl is decoupled into separate 
acceleration and braking rules Rla and Rib, after which their order is changed to 
Rib, R3, R4, Rla. The upshot of this is that there are no stopped vehicles in the 
system, thereby reducing the number of possible states for a cell by one. If v max = 1 
cell/time step, then the system can be fully described by cell occupancies. Apply- 
ing this SOMF theory to the STCA model, results in considerably underestimation 
of the flow in the system (even for the restricted case of t> max = 1 cell/time step) 
[140,72,54]. 

5.3.2 Improving the SOMF theory 

As mentioned in the previous section, setting v max = 1 cell/time step leads to an 
underestimation of the flow. However, when switching from a parallel update pro- 
cedure to a random sequential one, the resulting SOMF theory becomes exact ! 
It turns out that the reason for the underestimation, can be traced back to its ne- 
glecting of all correlations between cells (which are a consequence of the parallel 
update procedure). As explained in the beginning of section 3.2.4, using a parallel 
update excludes certain Garden of Eden states. However, the SOMF theory naively 
includes these paradisiacal states. As a solution, these GoE states can be eliminated, 
resulting in a paradisiacal mean-field theory (pMFT). In systems with higher max- 
imum speeds, more GoE states occur, making it difficult to derive a pMFT. Even 
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then, the theory still remains an approximation (albeit a better one) when using a 
parallel update procedure [141,72,54]. 

Taking into account short-range correlations, can be done by considering a car- 
oriented mean-field theory (COMF). Instead of dealing with cells and their occu- 
pancies, the COMF theory computes the probabilities P n (v) of finding a space gap 
of n cells for a vehicle driving with speed v [142]. In a sense, the COMF the- 
ory approximates the problem by neglecting the correlations between space gaps 
of successive vehicles [3]. As such, it gives qualitatively good approximations for 
p — > 0; in all other cases, the COMF theory starts to fail, because there are also 
correlations between the space gaps [18,3]. Note that the COMF theory has also 
been applied to the BJH-TCA and VDR-TCA models (see sections 3.3.2 and 3.3.3, 
respectively) [93]. 

Another approach to analytically solve the master equation, is to explicitly take 
into account the correlations between neighbouring cells, by considering clusters 
composed of n consecutive cells [72,54]. Such a site-oriented cluster-theoretic ap- 
proach proves to perform better than the COMF theory from the previous sec- 
tion [140]. The improvement of the approximation is even better when considering 
larger clusters; it is exact for n — > +oo [143,18,3]. 



6 Summary and outlook 

This report gave an elaborate and understandable review of traffic cellular automata 
(TCA) models, which are a class of computationally efficient microscopic traffic 
flow models. TCA models arise from the physics discipline of statistical mechan- 
ics, having the goal of reproducing the correct macroscopic behaviour based on a 
minimal description of microscopic interactions. 

We began with an overview of cellular automata (CA) models, their background 
and physical setup. Applying this technique to the modelling of traffic flows, we 
discretise a road into a number of small cells (a procedure called coarse graining), 
having a width of e.g., AX = 7.5 m. Time is also discretised into units of approxi- 
mately AT = 1 s. After introducing the mathematical notations, we showed how to 
perform measurements on a TCA model's lattice of cells, and how to convert these 
quantities into real-world units and vice versa. 

Subsequently, we gave an extensive account of the behavioural aspects of several 
TCA models encountered in literature. Already, several reviews of TCA models 
exist, but none of them consider all the models exclusively from the behavioural 
point of view. In this respect, our overview fills this void, as it focusses on the be- 
haviour of the TCA models, by means of time-space diagrams, (k,q) diagrams and 
the like, and histograms showing the distributions of vehicles' speeds, space, and 
time gaps. In the report, we have distinguished between single-cell and multi-cell 
models, whereby in the latter vehicles are allowed to span a number of consecutive 
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cells. We concluded with a concise overview of TCA models in a multi-lane set- 
ting, and some of the TCA models used to describe city traffic as a two-dimensional 
grid of cells, or as a road network with explicitly modelled intersections. The final 
part of the report illustrated some of the more common analytical approximations 
to single-cell TCA models. 

Considering the state-of-the-art in using TCA models, our analysis indicates that 
the field has evolved rapidly over the last decade. Starting from initial attempts 
based on rather crude models, the past few years have seen an increase in the com- 
putational complexity as well as the available computational power. More complex 
models are developed, of which we believe the brake-light TCA model of section 
4.2.2 is the most promising: it is able to faithfully reproduce the correct real-life em- 
pirical observations, and quite some work has been done at calibrating the model, 
see e.g., the recent work of Knospe et al. [19]. To conclude, we note an evolving 
trend of using these TCA models as the physical models underlying multi-agent 
systems, in part describing the behaviour of individual people in large-scale road 
networks [1]. 



A TCA+ Java™ software 

As already briefly mentioned in the paper, all simulations were performed by means 
of our Traffic Cellular Automata + software [102]. It was developed for the Java™ Vir- 
tual Machine (JVM), and can be downloaded 1 from: 

http : / / smtca . dyns . cx 

The software is also referenced on the Traffic Forum 2 (see section Links, subsec- 
tion Online Traffic Simulation or Visualization (Java Applets), item Java (Swing) 
application for several cellular automata models). 

In this appendix, we summarise our rudimentary TCA+ software. We start with an 
overview of its features, explain how to run the software, and conclude with some 
technical details with respect to the implementation of its code base. 

A. 1 Overview and features 

The TCA+ software package's goal is two-fold: on the one hand, it provides an 
intuitive didactical tool for getting acquainted with the concept of single-lane traffic 
cellular automata models. On the other hand, it provides a rich enough code base 
to perform hand-tailored simulation experiments, as well as giving insight into the 
details of programming TCA models. 

1 From May 2002 until June 2005, the software has been downloaded some 800 times, of 
which we suspect one third to be traffic coming from search engines' indexing robots. 

2 http://www.trafficforum.org 
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In a nutshell, our software considers one-dimensional traffic cellular automata with 
periodic boundary conditions, i.e., vehicles driving on a unidirectional circular 
road. Different sets of rules can be chosen, and for each set its parameters (e.g., 
stochastic noise) can be changed at run time. Both local and global measurements 
can be performed on the lattice by means of artificial loop detectors. A traffic light 
with cyclical red and green phases was also added, allowing to study elementary 
queueing behaviour. In the software, we have implemented the TCA models listed 
in Table A.l. 



TCA model 


Refer to section 


CA-184 


3.1.1 


DFI-TCA 


3.1.2 


STCA 


3.2.1 


STCA-CC 


3.2.2 


SFI-TCA 


3.2.3 


TASEP 


3.2.4 


ER-TCA 


3.2.5 


deterministic T 2 -TCA 


3.3.1 


stochastic T 2 -TCA 


3.3.1 


VDR-TCA 


3.3.3 


VDR-CC-TCA 


3.3.3 


TOCA 


3.3.4 


MC-STCA 


4.1 


HS-TCA 


4.2.1 


BL-TCA 


4.2.2 


KKW-TCA 


4.2.3 



Table A.l 

All TCA models implemented in our TCA+ software, accompanied by references to the 
respective sections in the paper where they are extensively discussed. 

In Fig. A.l, we show a screenshot of the main graphical user interface (GUI). As 
can be judged from the image, the TCA+'s GUI is rather huge, spanning approx- 
imately 1400x1200 pixels (scrollbars are automatically placed if it does not fit on 
the screen). It consists of several panels: 

• a scrolling time-space diagram containing vehicle trajectories and an animation 
of the road situation, 

• a panel containing some simulation statistics, 

• several simulator controls, 
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• and scrolling loop detector plots and plots of the (k,q), (k,v s ), and (q,v s ) dia- 
grams. 

In the following paragraphs, we describe each of these features in more detail. 
Note that there currently are two versions of the GUI: a standard version for all 
the single-cell TCA models, and a modified multi-cell TCA version with limited 
functionality (mainly for creating coloured tempo- spatial diagrams). 




Fig. A.l. A screenshot of the TCA+'s main graphical user interface (GUI) for single-cell 
TCA models. The GUI is rather huge, spanning approximately 1400x1200 pixels, con- 
sisting of several panels: a scrolling time-space diagram containing vehicle trajectories, 
an animation of the road situation, a panel containing some simulation statistics, several 
simulator controls, scrolling loop detector plots and plots of the (k,q), (k,v s ), and (q,v s ) 
diagrams. 

Vehicle animation 

Looking at the time-space diagram in the upper-left panel, we can discern the indi- 
vidual vehicle trajectories, as well as the typical backwards-travelling shock waves 
of congestion. In this scrolling diagram, the time axis goes from the left to the 
right, while the space axis goes from the bottom to the top (and is a one-to-one 
mapping of the consecutive cells on the ring road). Each pixel here corresponds 
to a unique cell of the simulator and each vehicle is coloured with a certain shade 
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of yellow (in order to easily distinguish between different neighbouring vehicles). 
There is also a setting available that allows stopped vehicles to be coloured red. 
In the upper-middle panel, the actual geometrical configuration of the ring road is 
depicted. This allows us to view the current physical situation on the road, i.e., the 
positions of all the vehicles. Each vehicle can be coloured with a certain shade of 
yellow (the same as in the time-space diagram). The current phase of the traffic 
light is also shown, as well as the positions of all the loop detectors: their positions 
are indicated by the small purple boxes alongside the road. The small green box in- 
dicates the position of the traffic light, with vehicles travelling in clockwise fashion. 

Simulation statistics 

In the upper-right panel, we can find the length of the ring road (expressed in the 
number of cells), the number of vehicles currently in the simulator, the global ve- 
hicle density, and the current time step. There is also a small panel that allows to 
quickly set the status of the traffic light to either red or green. 

Simulator controls and settings 

The middle-left panel contains buttons for starting, stopping (i.e., pausing), re- 
setting, and quitting the simulator. Several preferences can also be specified, i.e., 
whether or not to activate several panels containing the simulator's output. There 
is also the possibility to log the measurements from the loop detectors to a default 
file (called detector-values. data). And finally, the type of traffic cellular automaton 
(i.e., its rule set) can also be selected from a list, specified by radio control buttons. 

Note that there are several initial conditions possible for each density level: it is 
possible to start with a homogeneous state (all vehicles are spaced evenly), with a 
compact superjam of vehicles that are all stopped, or with a random initialisation 
(see also the introduction of section 3). 

If the simulation goes (visually) too fast, the cycle hold time can be increased, 
thereby freezing the simulation for a while between two consecutive time steps. 
Besides this, the ring road's global density and the vehicles' maximum speed can 
be specified. The sampling time for the artificial loop detectors can be adjusted (to 
increase or smooth out fluctuations). And finally, all probabilities can be adjusted 
between 0% and 100% in incremental steps of 1%. 

The red and green cycle times for the traffic light can be specified, such that the 
light can operate automatically, thereby inducing artificial queues at regular in- 
tervals. One can also control the traffic light manually (enabling the red or green 
phase) using the small upper-right panel; but if applied, the traffic-light controls 
override these manual settings. 

Plots of macroscopic measurements 

The software has the ability to extract both local and global macroscopic flow mea- 
surements from several uniformly road-side placed loop detectors which record 
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flows, densities, and space-mean speeds. 

The three large coloured regions in the middle panel represent the measured (and 
averaged) values of the local flows, local densities, and local space-mean speeds 
of the loop detectors. Pair-wise correlating these values, results in the plots of the 
(k,q), (k,v s ), and (q,v s ) diagrams in the lower-right panel. The coloured dots in- 
dicate locally obtained measurements, whereas the black dots represent globally 
obtained ones. 

Note the small button that allows to construct these diagrams: when it is pressed, 
the global density is incrementally increased from 0% to 100%, each time adding 
a single vehicle to the ring road. The simulation is then ran for a certain amount 
of time and the measurements from all the loop detectors are recorded. When all 
densities are processed (an indicator of the total time left is shown), the diagrams 
should be clearly visible in the loop detector plots in the lower-right panel. 

A.2 Running the software 

When visiting the website mentioned in the introduction of this appendix, there are 
two options for downloading the software. One is by downloading the compiled 
classes, whereas the other is to download the programme's source code. Once the 
compiled software has been downloaded, it is relatively easy to start the graphi- 
cal user interface. Considering the single-cell setup GUI, the software is ran by 
executing the following command: 

java -jar tea. jar 

Note that a Java™ Development Kit (JDK) (preferably Sun's 3 ) should be installed. 
Furthermore, due to a change in the threading of the Java™ Swing™ API, it appears 
that only JDK/JRE 1.3.1 is suitable ! 

A. 3 Technical implementation details 

It should be noted that the software is not implemented as an applet, but instead as 
a full Java™ application because it uses Swing™ components that are not standard 
supported by most browsers (at least not without installing a necessary plugin). The 
source itself logically consists of three different parts: 

• the TCA engine with different rule sets, 

• the graphical user interface, 

• and a whole range of predefined experiments. 

The geometrical configuration used in the single-cell TCA engine is a unidirec- 
tional ring road with a single lane. Vehicles are located in cells of AX = 7.5 m and 
can have speeds of to 5 cells/time step (corresponding to a maximum speed of 
135 km/h). One iteration in the simulation corresponds to a time step of AT = 1 s. 

3 http://java.sun.com 
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A number of artificial loop detectors are uniformly placed alongside the road, 
aggregating various macroscopic traffic measurements (i.e., flows, densities and 
space-mean speeds). In the GUI, global measurements on the entire lattice are per- 
formed according to the methodology explained in section 2.3.2, whereas local 
measurements are performed according to section 2.3.1. Note that for the TCA 
software itself, it is also possible to perform local measurements using a detector 
of unit length, according to the methodology explained in section 2.3.3. 



Besides the standard single-cell GUI and the limited multi-cell GUI, there also ex- 
ist some predefined experiments. These allow to create the (k,q), (k,v s ), and (q,v s ) 
diagrams, histograms of the vehicles' speeds, space gaps, and time gaps, as well as 
several order parameters (density correlations, nearest neighbours, and an inhomo- 
geneity measure that compares the locally recorded densities to the current global 
density). 



Inside the TCA+ software, several packages are available: 



• t ca . base containing the definitions of cells, global states, loop detectors, and 
the traffic cellular automaton's lattice, 

• tea. automata containing implementations of all the TCA models mentioned 
in section A.l, 

• tca.simulator containing the classes related to the single-cell and multi-cell 
GUIs, 

• tea . experiments . f undamentaldiagrams, 
tea . experiments . histograms, 

and tea . experiments . orderparameters containing setups for the pre- 
viously mentioned experiments. 
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B Glossary of terms 



B.l Acronyms and abbreviations 



ASEP asymmetric simple exclusion process 

BCA Burgers cellular automaton 

BJH Benjamin, Johnso, and Hui 

BJH-TCA Benjamin- Johnson-Hui traffic cellular 
automaton 

BL-TCA brake-light traffic cellular automaton 
BML Biham, Middleton, and Levine 

BML-TCA Biham-Middleton-Levine traffic cellular 

automaton 
CA cellular automaton 

CA- 1 84 Wolfram's cellular automaton rule 1 84 
ChSch-TCA Chowdhury-Schadschneider traffic 

cellular automaton 
CML coupled map lattice 

COMF car-oriented mean-field theory 

DFI-TCA deterministic Fukui-Ishibashi traffic 

cellular automaton 
DLC discretionary lane change 

ECA elementary cellular automaton 

ER-TCA Emmerich-Rank traffic cellular 

automaton 
GoE Garden of Eden state 

HS-TCA Helbing-Schreckenberg traffic cellular 

automaton 
JDK Java™ Development Kit 

KKW-TCA Kerner-Klenov-Wolf traffic cellular 

automaton 
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LGA lattice gas automaton 

LWR Lighthill, Whitham, and Richards 

MC-STCA multi-cell stochastic traffic cellular 

automaton 
MFT mean- field theory 

MLC mandatory lane change 

MPA matrix-product ansatz 

NaSch Nagel and Schreckenberg 
NCCA number conserving cellular automaton 
OVF optimal velocity function 

OVM optimal velocity model 

PCE passenger car equivalent 

PCU passenger car unit 

pMFT paradisiacal mean-field theory 
SFI-TCA stochastic Fukui-Ishibashi traffic 

cellular automaton 
SMS space-mean speed 

SOC self-organised criticality 

SOMF site-oriented mean-field theory 
SSEP symmetric simple exclusion process 

STCA stochastic traffic cellular automaton 

STCA-CC stochastic traffic cellular automaton 

with cruise control 
T 2 -TCA Takayasu-Takayasu traffic cellular 

automaton 

TASEP totally asymmetric simple exclusion 
process 

TCA traffic cellular automaton 

TMS time-mean speed 
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TOCA time- oriented traffic cellular 

automaton 

TRANSIMS Transportation ANalysis and SIMulation 
System 

UDM ultra-discretisation method 

VDR-TCA velocity-dependent randomisation traffic 
cellular automaton 
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B.2 List of symbols 



C(0) a CA's initial configuration 

C(t) a CA's global configuration at time step t 

5 a CA's local transition rule 

G a CA's global map 

G -1 a reversible CA's inverse global map 

K c the number of cells in one lane of a TCA's lattice 

C a CA's lattice (e.g., Z 2 ) 

Mi the (partially) ordered set of cells in the 

neighbourhood of the i th cell 
\J\f\ the number of cells in the neighbourhood 

of each cell 

Oc(t)\G- [ me backward orbit of the configuration C(t) 
under G -1 

®c(o)\G me f orwar d orbit of the initial configuration 

C(0) under G 
ai(t) the state of the i th cell at time step t 
E the set of all possible states a CA's cells can 

be in (e.g., Z 2 ) 
S £ the set of all possible global configurations 

of aCA 

T/^ the set of all possible configurations of a 

cell's neighbourhood 
| Tp M I the number of all possible rules of a CA 
Tc(o)\g the trajectory /orbit of the initial configuration 

C(0) under G 

a the entry rate of particles in the TASEP model 

ttj the anticipatory driving parameter of vehicle % 

a the acceleration capability of a vehicle in the 

KKW-TCA model 
(3 the exit rate of particles in the TASEP model 

b the deceleration capability of a vehicle in the 

KKW-TCA model 
bi (t) the state of the brake light of vehicle i at time 

t in the BL-TCA model 
5 the probability for a particle to move to the 

right in the TASEP model 
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A acCi the deterministic acceleration of vehicle % in 

the KKW-TCA model 
AT a TCA's temporal discretisation 
AV a TCA's speed discretisation 
AX a TCA's spatial discretisation 
D a parameter for the synchronisation distance 

in the KKW-TCA model 
D\ a parameter for the synchronisation distance 

in the KKW-TCA model 
Di the synchronisation distance of vehicle i in 

the KKW-TCA model 
r]i the stochastic acceleration of vehicle i in 

the KKW-TCA model 
7 the probability for a particle to move to the 

left in the TASEP model 
g s the average space gap 

g*. (t) the effective space gap of vehicle i at time t 

in the BL-TCA model 
fl'ssecurity a security constraint for the space gap in the 

BL-TCA model 
g t the median time gap 

~g ta the safe time gap in the TOCA model 
h the upper limit to the interaction horizon in 

the BL-TCA model 
£(t) a random number in [0, 1 [ drawn at time t 

from a uniform distribution 
kg the global density of a TCA's 

lattice 

ki the local density of a TCA's lattice 

Kc the number of cells in one lane of a TCA's 

lattice 
£ a TCA's lattice 

U the length of vehicle % 

I the average length of all vehicles on a TCA's 

lattice 

L c the number of lanes in a TCA's lattice 
Mg s ., Vi the gap-speed matrix of the ER-TCA model 



88 



p the slowdown probability in [0, 1] 

p the slow-to-start probability in [0, 1] 

p a the acceleration probability in [0, 1] in the 

KKW-TCA model 
p ai a parameter for the acceleration probability in 

the KKW-TCA model 
p a2 a parameter for the acceleration probability in 

the KKW-TCA model 
p acc the acceleration probability in [0, 1] in the 

TOCA model 
p the braking probability in [0, 1] in the 

BL-TCA model 

the deceleration probability in [0, 1] in the 

KKW-TCA model 
Pd the slowdown probability in [0, 1] in the 

BL-TCA model 
p dec the deceleration probability in [0, 1] in 

the TOCA model 
p s the slow-to-start probability in [0, 1] in 

the BJH-TCA model 
p t the slow-to-start probability in [0, 1] in 

the T 2 -TCA model 
P n (v) the probabilities of finding a space gap of 

n cells 

for a vehicle driving with speed v 
q g the global flow of a TCAs lattice 
qi the local flow of a TCAs lattice 
t q . the interaction horizon in the BL-TCA model 
v desi the desired speed of vehicle i in the 

KKW-TCA model 
v p a parameter for the acceleration probability 

in the KKW-TCA model 
v SB the space-mean speed in the free-flow regime 
v Sg the global space-mean speed of a TCAs lattice 
v Sl the local space-mean speed of a TCA's lattice 
x l - b the longitudinal position of vehicle «'s 

left-back neighbour 
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x- s the longitudinal position of vehicle i's 

left-front neighbour 

the longitudinal position of vehicle i's 

right-back neighbour 
x['* the longitudinal position of vehicle i's 

right-front neighbour 
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